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1 Introduction 



Light exotic nuclei are available as secondary beams at various radioactive 
beam facilities all over the world. These unstable nuclei are generally weakly 
bound with few, if any, bound excited states. They have been studied ex- 
tensively in recent years by electromagnetic excitation with the help of the 
Coulomb breakup method [1,2,3,4,5]. For low orbital angular momenta of the 
lowest bound valence nucleon an extended diffuse density distribution, a halo, 
develops resulting in a large size of the nucleus [6,7,8,9,10]. Simultaneously, 
electromagnetic transitions to the continuum with large strength are observed 
at low energies. Properties of stable nuclei have successfully been investigated 
by electromagnetic excitation in photonuclear reactions as well as in heavy 
ion collisions for a long time. Their electromagnetic strength functions are 
dominated by the giant resonances high in the continuum. 

Nuclei close to the neutron and proton drip lines often exhibit a pronounced 
nucleon+core structure that is well described by single-particle models with 
appropriately chosen potentials. In a microscopic shell- model study [11] strong 
low-lying dipole strength in neutron-rich ^^Be and proton-rich ^^O was ob- 
served to originate from loosely-bound extended single-particle wave functions. 
These extended wave functions were obtained by adjusting the potential depth 
in order to reproduce the empirical binding energies. In order to study low- 
lying strength theoretically, such a feature has to be added to the ab-initio 
microscopic approaches. A comparison of strength functions deduced from ex- 
periment with theoretical predictions from single-particle models is used to 
extract spectroscopic factors or asymptotic normalization coefficients (ANC). 
They can be compared to more elaborated nuclear models. However, the im- 
portance of effects from the interaction between the nuclei in the final state 
has to be assessed in order to obtain reliable information from experimental 
data. Both the shape and the absolute magnitude of the strength function can 
be affected. This will also have consequences for the application of sum rules 
that relate the total excitation strength to the properties of the ground state. 

The relevant matrix elements for the electromagnetic transition to the con- 
tinuum at low energies are esssentially determined by the asymptotics of the 
bound-state and continuum wave functions of halo nuclei. This allows to study 
systematically the effects of the final-state interaction (FSI) without the ne- 
cessity to introduce sophisticated nuclear structure models. Magnetic contri- 
butions to the continuum transitions are usually much weaker than electric 
transitions, except for the deuteron breakup at low energies [12] and for the 
excitation of resonances. We will limit ourselves to the discussion of direct elec- 
tric transitions to the continuum, with emphasis on El transitions, but the 
search for low-lying Ml strength and its theoretical description [13] remains 
an interesting challenge for future studies. 
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The nuclear interaction Vbc between the nucleon b and the core c in the elec- 
tromagnetic breakup of an exotic nucleus a is responsible for the binding of 
the nucleus a. However, it also affects the structure of the continuum (usually 
in partial waves with different I values) even if there are no resonances ob- 
served at low excitation energies. In the experimental analysis of neutron+core 
breakup reactions it is often neglected and a plane wave is assumed in the final 
state of the b + c system. The interaction between b and c also appears in the 
final state of the photo- dissociation reaction 0(7, c)b or in the initial state of 
the radiative capture reaction b{c, 7)0. Thus, the interaction Vbc can affect the 
energy dependence of the astrophysical S factor of the radiative capture reac- 
tion that is used to extrapolate experimental data to zero energy. Similarly, 
the strength of the interaction enters into the calculation of the zero-energy 
S factor from asymptotic normalization coefficients (ANCs) that are deter- 
mined experimentally from transfer reactions in the ANC method, see, e.g., 
[14,15,16,17]. 

Experimentally observed excitation functions of exotic nuclei show an approx- 
imately universal shape when plotted as a function of appropriately scaled 
variables [18,19,20,21]. They are dominated by direct transitions to the con- 
tinuum with only few resonances and simple scaling laws apply. The nuclear 
structure in the initial and final state depends only on a limited number of rel- 
evant quantities that contain all the structural information that is accessible 
at low energies. Details of the nuclear interaction are not resolved at this low 
energy scale. The interaction between the fragments leads to a change of the 
transition strength when compared to the case without nuclear interaction. 
The actual nuclear potential is often not well constrained since extrapolations 
of the corresponding systematic optical potentials, e.g. [22], from nuclei in the 
valley of stability to unstable nuclei are qiiestionable. The effect of the contin- 
uum interaction was studied before only in selected cases. E.g., it was found 
that the s-wave ground state to p-wave continuum El transition in the case 
of ^^Be is much less affected by the potential in the final state than the p — > s 
transition in [18,19]. 

At low energies the effect of the nuclear potential is conveniently described 
by the effective-range expansion [23,24], with the scattering length and the 
effective range as the main parameters. An effective-range approach for the 
FSI in electromagnetic excitations was introduced in [25] and applied to the 
breakup of ^^Be. Recently, the same method was applied to the description of 
electromagnetic dipole strength in ^^O [26]. Here, we will discuss this approach 
in much more detail. A systematic study for various transitions will shed 
some additional light on the sensitivity to the interaction in the continuum. 
We want to expose the dependence on the binding energy of the nucleon 
and on the angular momentum quantum numbers. Our approach extends the 
familiar textbook case of the deuteron [13], that can be considered as the 
prime example of a halo nucleus, to arbitrary nucleon-|-core systems. 
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Our effective-range approach is closely related to effective field theories that 
are nowadays used for the description of the nucleon-nucleon system and halo 
nuclei [27,28,29]. (The deuteron, the only bound state in the nucleon-nucleon 
system, can be considered as a good example of a halo nucleus). The charac- 
teristic low-energy parameters are linked to QCD in systematic expansions. 
Similar methods are also used in the study of exotic atoms {7i~A, tt+tt", 7i~p, 
. . . ) in terms of effective-range parameters in Ref. [30]. The close relation of 
effective field theory to the effective-range approach for hadronic atoms was 
discussed in Ref. [31]. In our approach these constants are treated as free 
parameters. It is not our aim to relate them to the underlying microscopic 
description. Aspects of the many-body physics are summarized, e.g., in terms 
of spectroscopic factors or asymptotic normalization coefficients. 

We also investigate in detail a specific model, the square well potential. It has 
great merits: it can be solved analytically, it shows the main characteristic 
features and it leads to rather simple and transparent formulas where some of 
them seem to be new. These explicit results can be compared to our general 
results for low energies (effective-range approach) and also to Woods-Saxon 
models. 

This paper is organized as follows. A nucleon-|-core potential model for halo 
nuclei is introduced in section 2. The relevant scaling parameters are defined 
in subsection 2.1 and scaling laws for the the probabilities to find the nu- 
cleon inside the range of the nuclear potential are discussed for bound and 
scattering states in subsection 2.2. The scahng of the root-mean-square ra- 
dius serves as another indication for the halo nature of the bound state. The 
reduced transition probabilities for the breakup of a nucleon-|-core nucleus 
arc calculated in subsection 2.3 for electric transitions with multipolarity A. 
They only depend on the asymptotic normalization of the bound state wave 
function and on radial integrals with the asymptotic wave functions since the 
radial integrals are dominated by the contribution from outside the nuclear 
potential. In the square-well model explicit expressions for the ratio of the 
interior to the exterior contributions are derived. An alternative calculation of 
El transition integrals with the help of a commutator relation is presented in 
subsection 2.4. The strength functions are related to cross sections of photo- 
nuclear reactions in subsection 2.5 and the high-energy behaviour is discussed. 
The effect of the nucleon-core potential is considered in subsection 2.6. The 
effective-range expansion allows to parametrize the effects of the continuum 
interaction in a suitable way to study systematically its influence on the tran- 
sition strength. In Section 3 the reduced transition probabilities are expressed 
in terms of characteristic shape functions that depend on certain reduced in- 
tegrals. The dependence of the shape functions on the scaling parameters is 
studied in various limits. Systematic variations of the shape functions are dis- 
cussed for neutron-|-core systems in subsections 3.1 and 3.2 without and with 
FSI, respectively, where analytical results are obtained. Proton-|-core systems 
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are treated numerically in subsection 3.3. The relation of the total excitation 
strength to the properties of the bound state with the help of sum rules is con- 
sidered in section 4. The findings of the model with asymptotic wave functions 
are corroborated in more realistic calculations using wave function generated 
from Woods-Saxon potentials in section 5. The implication of the nucleon-core 
interaction in the continuum state on the ANC method is discussed in section 
6. We close with a summary and conclusions. The appendix contains detailed 
derivations and explicit expressions of our analytic calculations. 



2 Nucleon-|-core model 



Exotic nuclei close to the driplines often exhibit a pronounced structure with 
a nucleon b (proton or neutron) weakly bound to a core c. They are often well 
described by simple single-particle models that are able to explain the ba- 
sic features of low-energy excitations. For small separation energies Sb of the 
nucleon and low orbital angular momenta I the exotic nucleus a develops a 
proton or neutron halo where the nucleon wave function extends to large radii 
and there is a large probability of finding the nucleon outside the classically 
allowed region of the potential Vbc- Matrix elements for electromagnetic tran- 
sitions only depend on a small number of characteristic scaling parameters. 



2. 1 Scaling parameters of halo nuclei 



The main scale is set by the nucleon separation energy Si, {b = n,p) that is 
related to an inverse decay length 

V^Sb 



of the bound state with the reduced mass fi = mbmc/{mb + nic). It becomes 
very small for typical halo nuclei with small binding energies. The scattering 
state is characterized by the momentum hk that is related to the relative 
energy in the continuum E by 



k = ^. (2) 



A third parameter is the size of the nuclear core given by the range of the core- 
nucleon potential R. These three parameters allow to define the dimensionless 
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quantities 



7 = qR and k = kR (3) 



and the ratio 

K k , ^ , , 

1 q 



independent of R. It can be considered as a definition of halo nuclei that 
the parameter 7 is small. This means that the extension of the wave function, 
characterized by is much larger than R. The parameter 7 can be used as a 
convenient expansion parameter in a systematic approach to calculate matrix 
elements. On the other hand, the relevant range of the parameter x extends 
from zero to a value of the order of one. For larger values of x other degrees of 
freedom, like core excitations, will come into play and tend to invalidate the 
simple model. 

For proton+core systems an additional scale enters the problem set by the 
Gamov energy 



or the nuclear Bohr radius 



with charge numbers Zi, and of the nucleon h and the core c. They serve to 
define the Sommerfeld parameters 




m^-^ = and 77/ = = (7) 

qaN V kttN 



of the bound (i) and the scattering (/) state with the relation 

xrjf = rji . (8) 

The nuclear Bohr radius (6) is related to another important parameter: half 
the distance of closest approach in a head on collision with energy E. It is 
given by = l/(/c^ajv). Typical features of a halo system essentially depend 
on the three independent parameters j, x, and rji. The characteristic scaling 
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Table 1 

Characteristic scaling parameters for electromagnetic strength in single-particle halo 
nuclei. 



origin energy scale dimensionless parameter 



bound state bf, = — — 

Z/jL 

one-nucleon separation energy 


7 = qR 




fi^k"^ K k 
scattering state E = — — x = — = — = \ 

2fi 7 q \ 

nucleon-corc relative energy with k = kh 


I 


Coulomb field Eq = | 


(Zi,Z,e^\'lx 

[ h ) 2 2„a% "^'-y 
Gamov energy 


[E^ _ 1 
Sb qaN 


= xr)f 



parameters are summarized in Table 1. For a neutron+core system one obvi- 
ously has r)i = 0. In case of a proton-|-core system with a large charge number 
of the core and/or small binding energy Sf, the parameter rji can become quite 
large. 



2.2 Probabilities in nucleon+core systems 



Common to all nucleon-|-core nuclei with small binding energy is the large 
probability in the bound state of finding the nucleon outside the range of 
the nuclear potential. A similar observation is made for the scattering state. 
For halo nuclei the miclcon does not penetrate strongly into the range of the 
nuclear potential that describes the halo nature of the bound state. 

In a simple neutron-|-core model assuming a square-well potential of radius R 
the probability Pni of finding the neutron with principal quantum number n (= 
number of nodes of the radial wave function including the node at r = 0) and 
orbital angular momentum / inside the range of the potential can be calculated 
analytically (see Appendix A). It essentially depends on the parameter 7. One 
finds that the typical halo structure appears only for low I, i.e. s, p waves, and 
small neutron separation energies S^, i.e. small 7. Figure 1 clearly shows the 
increase of the probability to find the neutron outside the range of the potential 
with decreasing 7. For larger values of / the centrifugal barrier hinders the 
occurence of a halo structure and the penetration into the classically forbidden 
region is reduced. A larger number of nodes in the wave function increases the 
halo effect again. This effect is most pronounced for s waves. In the extreme 
halo limit 7^0, however, there is no dependence on n any more and the 
probability approaches finite values of 0, 1/3, and 3/5 for s, p, and d waves, 
respectively. In case of a more realistic Woods-Saxon shape of the nuclear 
potential the probability is even smaller as in a square-well potential with the 
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Fig. 1. Probability P^i of finding a neutron with separation energy 5„ inside the 
radius i? of a square- well potential as a function of the parameter 7^ = ijjiSnB? /f? 
for orbital angular momenta i = 0, 1, 2 and principal quantum numbers n = 1 (thick 
line) and n = 2 (thin line). 

same radius [32]. This is easily understood since the depth of the potential is 
reduced inside the radius and increased outside the radius. As a consequence 
the probabihty is shifted to larger radii. 

A different way to characterize the halo effect for small orbital angular mo- 
menta is the dependence of the root-mean-square (rms) radius \J~(t^i on the 
parameter 7. Explicit expressions of the rms radius for / = 0, 1, and 2 in 
the square- well model were given in Ref. [33], see also [32]. However, simpler 
expressions for arbitrary values of I can be obtained (see Appendix A). The 
scaling behaviour of the rms radius is well known, see, e.g., [6,7,8,9,34,35,36] 
and Figure 2. For / = and / — 1 one finds the scaling laws 



^ if / = 1 

for 7 — > 0. For / > 2 the rms radius approaches a finite value. The divergence 
of (r^); for / = and 1 is the typical sign of the halo nature. 

Large values for I and small parameters 7 prevent the neutron in a continuum 

state to enter into the nuclear interior as long as there is no resonance. In 
figure 3 the corresponding differential probability dP^i/dx (see Appendix A) 
as a function of the ratio x = k/j is shown for various parameters 7 and 
The explicit expression (A. 37) shows a scaling proportional to 7 and inversely 
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Fig. 2. Root-mean-square radius {r'^)^^i of a neutron with separation energy 5„ in 
a square-well potential of radius i? as a function of the parameter 7^ = 2iiSn,Fl? /h^ 
for orbital angular momenta Z = 0, 1, 2 and principal quantum numbers n = 1 (thick 
line) and n = 2 (thin line). 
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Fig. 3. Differential probability dPni/dx as a function of a; = k/j of finding a neu- 
tron with orbital angular momenta Z = 0, 1,2 inside the radius i? of a squarc-wcU 
potential that binds a neutron with separation energy Sn, orbital angular momen- 
tum I and principal quantum numbers n = 1 (thick lines) and n = 2 (thin lines) 
for various values of the parameter 7^ = 2/xS'„i?^/^^ = 0.1 (solid line), 0.5 (dashed 
line), 1.0 (dotted line), 2.0 (dot-dashed line). 
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proportional to the penetrability. For larger x where resonances in the scatter- 
ing occur the probabiHty exhibits clear maxima. However, the neutron does 
not penetrate deeply into the nuclear interior for small x, small 7 and large /. 
Again, a larger number of nodes in the bound-state wave function that deter- 
mines the potential depth also reduces the probability to find the neutron in 
the potential well for continuum states. 

In the case of proton+core systems the probabilities of finding the nucleon in- 
side and outside the range of the potential, respectively, cannot be calculated 
analytically in the model with a square-well potential. Obviously, some quali- 
tative changes are expected. Due to the additional Coulomb barrier the proton 
penetrates less deeply into the classically forbidden region and the probability 
for the proton to be inside the nuclear radius is enlarged for bound states. 
On the other hand, a proton in a continuum state will be found less proba- 
ble inside the range of the nuclear potential as compared to a corresponding 
neutron with the same energy. This is described by the Coulomb penetration 
factor. 

Since the nucleon in halo systems can be found predominantly outside the 
range of the nuclear potential one can expect that the relevant transition ma- 
trix elements mainly depend on the asymptotics of the bound and scattering 
wave functions. This is especially true for electric multipole transitions that 
contain an additional dependence enhancing contributions from large radii. 

2.3 Reduced transition probabilities and radial integrals 

The reduced transition probability dB{7rX) /dE of multipolarity ttA (tt = E, M; 
A = 1,2,...) for the electromagnetic breakup of the nucleus a into b + c with 
relative energy E is the basic quantity of our study. It contains the information 
on the nuclear structure in the initial ground state and the interaction in the 
final continuum state. This strength function determines the response of the 
system to a photon and enters the expressions for the corresponding cross 
sections. It can be extracted from experimental data in order to compare 
directly to predictions of nuclear models. 

Assuming the nucleon-|-core picture for the halo nucleus the reduced transition 
probabilities are easily expressed in terms of certain radial integrals with the 
wave functions in the initial (bound) and final (scattering) states denoted by 
i and /, respectively, in the following. In general, the spin s = 1/2 of the 
nucleon b couples with the orbital angular momentum k/f to the total angular 
momentum j'j/j. The total angular momentum Jj/j of the system a — b + c is 
obtained by coupling j^/f with the spin of the core jc- 

Usually there are various combinations of ji/f and jc that are possible to 
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contribute to a given value of Jj/j. The reduced transition probability for a 
specific electromagnetic transition ttA to a final state with momentum hk in 
the continuum is given by 



dB 

— (ttA, JiS kJfs) = 



(10) 



2J^ 
2 Ji + 1 ^ 



{kJfjflfsjcWMinX) 1 1 Jijiksjc) 

jihjc 



fj,k 



depending on reduced multipole matrix elements. In the following we will only 
consider electric excitations (tt = E) with multipole operator 



:iii 



that dominate the continuum breakup of exotic nuclei. The effective charge 
number is given by 

Z^£> = ZJ-^)\Z.U^)\ (12) 



For proton+core systems the effective charge numbers for El and E2 transi- 
tions are of comparable magnitude and, generally, one has to consider both 
contributions in the cross sections for Coulomb breakup, photo dissociation 
or radiative capture. In the case of a neutron+core system, the E2 effective 
charge number is suppressed by a factor 1/A as compared to El since = 0. 
El transitions dominate the low-lying electromagnetic strength and the E2 
contribution can be neglected. Neverthess we include the E2 case here for com- 
pleteness. The strong reduction of contributions from higher multipolarities 
in the neutron-|-core case was noticed before, e.g. in Refs. [21,37]. 

In the single-particle model the wave functions of the initial and final state 
are given by 



^iir)^{r\JijiliSjc 



(13) 



and 



^fir) = {r\kJfjflfsj,) (14) 
An 
kr 



^ E (j/^/Jcme|J/M^)4^.^,^(r)z^/F,;^/^)J^jj;^(r)0,.„.., 



m/mc 
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respectively, with the radial wave functions fJ^j^l^i'>^) and g'j^j^i^ir) and with 
the spinor spherical harmonics 



= ^{Ims ms\j m)Yi^{f)xs 



(15) 



minis 



The wave function of the core is denoted by 4)j^mc - The reduced matrix element 
in (10) can be expressed as 



{kJfjflfsjc I \M{EX) 1 1 Jijiksjc 



k 



-^^^(A^^c) (-^)'^/gi:^(Ajc) 



with the angular momentum coupling coefficient 

= (-l)^+^*+'/+^(-l)^=+-^+^/+^(/, A 0\lf 0) 
^J2ji + 1^2/^ + 1^2 Ji + Ipjf + 1 

/ 2A + 1 
V 47r 

and the radial integral 



^ ji 

< 


> < 


ji jc Ji 

> 









'■Jd. 



t'{Xjc}-Jdrg^;^^,^iryf^^,^{r) 



(16) 



(17) 



(18) 



that contains the radial wave functions of the bound state fj^j^^ir) and the 
continuum state g'j^j^i^ir), respectively. 



The asymptotic of the radial wave functions for the ground state 



Ji 1 



(19) 



is determined by the asymptotic normalization coefficient Cj^^.f. of the true 
many-body wave function and a Whittaker function VF^_^.^i.+i/2 [38]. The 
bound state is characterized by the parameters q, rji, and the orbital an- 
gular momentum Zj. In the pure single-particle model for the nucleon-|-core 
system the corresponding ANC Cj':jj_{sp) is determined by the normaliza- 
tion of the wave function /j^j-;. ('")• Because the radial wave function for small 
r depends on the nuclear potential of the single-particle model there is a 
model-dependence of Cj^^-.;.(sp). Similarly, the spectroscopic factor Sj:j^i, for 
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the nucleon+core configuration depends on the single-particle model that is 
used for the calculation. In contrast, the actual ANC of the true many-body 
wave function 



(J3c _ (Jjc / \ ^Jc 



(20) 



is a model-independent quantity. It is directly inferred from transfer reactions 
for example. In the following we always assume a spectroscopic factor of one, 
i.e. the single-particle ANC and the true ANC are identical. 

For the scattering state we have the asymptotic form 



with regular and irregular Coulomb wave functions Fi^. and G/^ [38], respec- 
tively, and Coulomb phase shifts aij^ that depend on the Sommerfeld parameter 
Vf = Vi/^- Effects of the nuclear interaction in the final-state are contained in 
the nuclear phase shifts ^^jjjjij- 

In general the radial integral (18) has to be calculated with the relevant bound 
and scattering wave functions with the correct asymptotics (19) and (21), re- 
spectively. They are obtained by solving the Schrodinger equation for a given 
potential, e.g. of Woods-Saxon form V{r) — —Vq/{1 + exp[(r — i?o)/a]) with 
depth Vo, radius Rq and diffuseness parameter a. In the following calculations 
we will only consider nuclear central potentials and neglect contributions from 
the spin-orbit interaction to simplify the discussion. In case of proton+core 
systems the contribution of the Coulomb potential can be assumed to be that 
of a homogeneously charged sphere of the same radius Rq as the nuclear poten- 
tial. The potential depth Vq has to be adjusted to give the correct separation 
energy Si, of the nuclcon in the ground state of a. However, in the continuum 
state the parameter Vq can be varied freely to investigate the dependence of 
the transition strength on the final-state interaction. 

A remarkably good approximation for the nuclear interaction is the square- 
well potential and many results can be derived analytically. They already 
show the main features for the scaling laws of the matrix elements. In this 
case the contributions to the integral (18) from the interior and exterior part 
can be calculated independently to estimate their relevance in the transition 
matrix element. In Appendix A. 4 explicit expressions for the dipole integrals 




(21) 
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Fig. 4. Modulus of the ratio -Rj^(l) of the interior and exterior radial integrals for 

dipole transitions s ^ p (top) and p ^ s (bottom) as a functions of x"^ = E/Sn for 
various values of the parameter 7^ = 2fj,SnR'^ /h? ■ Thick (thin) lines correspond to 
the pricipal quantum number n = 1 (n = 2). 

(suppressing irrelevant quantum numbers) 

R 00 

ll^il <) = j dTg;^{T)rfi^{T) and >) = j dr g^^^ry fi^r) (22) 



R 



in the ncutron+core case for s p and p ^ s transitions are derived. They 
only depend on the value of the radial wave functions and their logarithmic 
derivative at the radius R of the square well. The ratio 



' ^(1>) 



(23) 
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as a function of for s — > p and p ^ s transitions, respectively, is displayed 
in figure 4 for various values of 7^ assuming the same depth of the square-well 
potential for the bound and scattering state. For halo nuclei with small 7^ 
the total radial integral is clearly dominated by the exterior integral even at 
values — E / Sn ^ I- I'a the case of a resonance in the scattering wave (e.g. 
for the s ^ p transition with 7^ = 1) the ratio b''/. (1) shows a distinct peak 
since the scattering wave function penetrates into the core of the nucleus. In 
the limit a; — > one finds the scaling laws 



depending on the principal quantum number n = 1,2,... (see appendix A. 4). 
For larger values of n the ratios become smaller but they remain finite for a 
given 7 in the limit a; — * 0. As can be seen from figure 4 these scaling laws are 
also well satisfied for larger values of x'^ as long as there are no resonances or 
accidental zeros in the interior integral. For halo nuclei one clearly sees that 
the total radial integral is well approximated by the exterior part with the 
asymptotic form of the wave functions that is independent of the details of 
the potential. This will be even more true for higher multipolarities due to the 
factor in the radial integral. Only the parameters 7, k, the ANC C;. and 
the phase shift are really relevant. 



2.4 Dipole integrals and commutator relations 



In the case of electric dipole transitions the relevant reduced matrix element 
(16) can also be calculated with the help of the commutator relation [39,40] 

{^f\ [H, [H,r\] 1$,) = {Ef - Eif{^f\i^^i) . (25) 
Introducing the scaling parameters 7 and x, we find 



- ^ir^7?^(^/l [H,r\] m (26) 



for the dipole matrix element. This equation directly shows the occurence of 
a pole at x"^ — —1, i.e., E — —Si,. This derivation is much more transparent 
than the corresponding discussion in [41] for El radiative capture reactions. 
For halo systems with small nucleon separation energies one immediately 
sees that a series expansion of the dipole matrix element (26) in the parameter 
X (or the energy E) is only of limited value because of the small radius of 
convergence Xconv — 1 (-E'conv = Sb). It would be more advantageous to expand 
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the matrix element with the double commutator [H, [H, r\] in a power series 
in X. 

The double commutator can be calculated assuming difTerent forms of the 
potential V in the Hamiltonian H ~ p^/(2/x) + V oi the system. The simplest 
case is a central potential V{r) that commutes with r, i.e. [V,r] = 0. With 
[H, [H, r\] — h^VV{r)/iJi the radial integral (18) can be expressed as 



In the neutron+core case with a square-well potential 

V{r) = -Voe{R - r) (28) 



we find that the radial integral is given by 

^a^^ (l^e) = j;^^,9%,,, im., iR) (29) 



It depends only on the values of the initial and final radial wave function at 
the radius R and the depth of the potential in the dimensionless parameter 

V - (30) 



that characterizes the strength of the interaction. It is not surprising that this 
result can also be derived by evaluating the radial integral (18) directly (see 
appendix A. 4). For more realistic potentials the derivative dV/ (dr) also peaks 
close to the nuclear radius and the radial integral is mainly sensitive to the 
bound and scattering wave functions near the nuclear surface. 

For finite values of R the radial integral (29) shows a x^^ dependence for large 
energies independent of the orbital angular momenta in the initial and final 
state. (The wave function gi. is independent of x and fi^ shows an oscillatory 
behaviour for large x.) This E'"^ dependence is sufficient for the convergence 
of the non energy-weighted and energy-weighted sum rules as discussed in 
section 4. 

In the proton-|-core case the potential is well approximated by 

Z 

V{r) = -Vo 9{R -r) + ^— e{r - R) (31) 
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and the radial integral becomes 



J ■ I 2i? 

^Jdik (lie) = ^^2 _^ ^2)2 (^2) 



X 

with the integral 

dr 

R 



that converges much more rapidly than the usual El transition integral. This 
contribution only depends on the parameters j, rji, and x. The ratio 

> iv + 2^r„)gt^iR)f,,iR) ^ > 



of the two contributions in (32) is depicted in Fig. 5 as a function of x'^ for 
typical values of 7^ and rji for both s ^ p and p ^ s transitions. It is clearly 
seen that the ratio Fi^{l) of the two contributions to the total radial integral 
increases with increasing rji and decreasing 7. This is expected since Coulomb 
effects become stronger for larger rji and the bound state wave function at 
larger radii is less steep for smaller 7. The ratio (34) decreases for larger 
values of x'^ with a larger effect for large values of 7. Generally, transitions 
p ^ s are more strongly affected by the correction than transitions s ^ p, 
however, for large 7 (less halo effect) the difference becomes smaller. 

Ususally, one can expect that the Hamiltonian H — p^(2/x) +V contains more 
general potentials that do not commute with r [42], e.g., 

V = Vc{r) + Vis{r)i- s (35) 



with a central potential Vc(r) and a spin-orbit contribution Vis{r). In this case 
the double commutator in (25) can be calculated explictly and a rather compli- 
cated radial integral is obtained. Alternatively, we can introduce Hj = PjHPf 
and Hi — PiHPi with projection operators Pf and Pf on the final and initial 
state, respectively. Then the commutator relation (25) can be generalized to 



{Ef - E,f{^f\T\^,) = ($^1 ([E, [E,r1] + [S, {A,f}] (36) 

+{A,[S,r-]} + {A,{A,f}})|$.) 
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2 



Fig. 5. Ratio of Eq. (34) as a function oi = E / Sp for various values of 

7^ = 2fiSpR'^ /h!^ and r]i for the bound state without a node (principal quantum 
number n = 1). Thick (thin) lines correspond to transitions s ^ p {p ^ s). 

where E = {Hf + Hi)/2, A = {Hf — Hi)/2 and {.,.} denotes the anticom- 
mutator. For the potential (35) we find — p^/(2/x) + Vfe(r) with k — i, f 
and 

Vk{r) = K(r) + UkUk + 1) - kik + 1) - s{s + 1)] (37) 



where jk, h, and s are the total angular momentum, the orbital angular mo- 
mentum, and the spin of the nucleon, respectively. The central potentials 
Vfe(r) with [Vfc(r),r] = are now different in the bound and scattering states, 
however, the various contributions in (36) can be calculated easily. Assum- 
ing a square- well potential for both Vc{R) and Vis{r) with the same radius 
but different depths the radial integral corresponding to (29) becomes a more 
complicated expression that contains also derivatives of the wave functions at 
the radius R. Explicit expresssions are obtained from adding the interior and 
exterior integrals as derived in appendix A. 4. Also other forms of the poten- 
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tials, e.g., a surface-peaked spin-orbit potential Vis = const, x 5{r — R), can 
be considered. However, we omit the details of the calculation. 



2.5 Cross sections 



The reduced transition probability determines the cross sections for photo- 
nuclear reactions. The photo dissociation cross section 



...(a + ,^6 + c) = — pA^y ^ (38) 



is proportional to dB{TTX)/dE where = E + Sh is the sum of the binding 
energy 5*;, > of the nucleus a with respect to the breakup into b + c and 
the relative energy E in the final state. The photo absorption cross section 
(38) also enters the cross section for the electromagnetic dissociation reaction 
a + X^b + c + X during the scattering of an exotic nucleus a on a target 
nucleus X. A first-order calculation gives 

d^a 1 ^ , , , driT^x 



E^-A(a + 7-& + c):^ (39) 



dO^axdE E^ dO^aX 



with virtual photon numbers driT^x/dVlaX that can be calculated in the semi- 
classical theory or in prior-form DWBA in the quantal approach. The fac- 
torization of the cross section (39) into contributions that are related to the 
nuclear structure of the exotic nucleus ct^a and to the excitation mechanism 
dn.j^x/ dftax is no longer valid if higher-order effects from the target-fragment 
interactions are significant. Prom the photo dissociation reaction it is also pos- 
sible to extract information on the radiative capture reaction 6 + c ^ a + 7 
that for low energies is relevant for nuclear astrophysics. The capture cross 
section 

a,.(6 + c ^ a + T) = p j'f^pj'V 1) + 7 ^ 6 + c) (40) 

is obtained by applying the theorem of detailed balance. The photo absorption 
cross section (38) is multiplied with a factor that contains the spins Jj of the 
particles i — a,b,c, the photon momentum hk^ — E^/c, and the final state 
momentum Hk. 

For electric dipole transitions it is easy to find the high-energy behaviour of the 
various cross sections. From equation (27) the dependence I^j (1) oc x'"^ cx E""^ 
of the radial transition integral is extracted. With (10) and (16) one finds 
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dB{El)/dE oc x-i|7;/(l)P oc x'^ or dB{El)/dE oc E'^l'^. This corresponds 
to a dependence o"si(a+7 &+c) oc E~'^/'^ and f7£;i(fc+c — > a+7) oc for 
the photo dissociation and radiative capture cross sections at high energiess 
E, respectively. We note that in the case of the photodissociation of atoms, 
with the 1/r shape of the Coulomb potential, one finds the same -E""^/^ law, 
see [43]. 



2.6 Effective-range expansion 



The nucleon-core interaction that is responsible for the binding of the system 
also generates structures in the continuum. At low energies in the continuum 
the phase shifts are insensitive to the details of the nuclear potential as long 
as no resonance appears. They are determined by only a few parameters that 
are given by the effective-range approximation. E.g., in the case of charged- 
particle scattering the low-energy phase shift 5i in the partial wave with orbital 
angular momentum / is determined by the scattering length ai and the effective 
range r; in the expansion [44] 



21+1 



277/1(77) 

cot 5; + 



-± + ^k' + ... (41) 



with the function 



/i(77) = l [^/>(1 +ir/) +^/;(l -^77)] -lnr7 (42) 

00 2 

+ 77"^) 

that depends on the Sommerfeld parameter 77 in the argument of the Digamma 
function tjj and 7 = 0.5772156649 ... is Euler's constant. The constants Cf for 
increasing I are obtained by means of a recursion relation 

Cf(,) = Ct:(,)(l4) with = (43) 



Note that ai and ri for Z > do not have the dimension of a length. The 
effective-range expansion can be used to express the phase shifts directly as 
functions of the momentum Hk. 

Taking only the scattering length ai and the effective range ri into account one 
does not necessarily find a good approximation of the phase shift 5i over a wide 
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range of energies in the continuum. However, the effective-range expansion 
motivates to introduce an energy-dependent function bi via the equation 



cot di + 



2rih{ri) 



(44) 



with the dimensionless parameter x — k/q where q is constant for a given 
nucleus. In general, the function bi depends on the momentum hk in the final 
state. It fully describes the effects of the interaction in the scattering state at 
all energies. The quantity bi varies slowly with the parameter x = k/q as long 
as no resonance is approached. For small k one can identify bi with a reduced 
scattering length since the usual scattering length is obtained in the limit 



21+1 



ai — lim — 
x^o\q 



(45) 



and for Si — we obviously have 6^ = 0. 



3 Reduced radial integrals and shape functions 



Considering the asymptotics of the radial wave functions one can define di- 
mensionless reduced radial integrals Tj^.^iJ i^jc) by the relation 



exp 



(46) 



At low energies the main contribution to the radial integral arises from radii 
r larger than the radius of the nucleus. This is especially true for halo nuclei 
where the probability of finding the nucleon inside the range of the nuclear po- 
tential is small (see subsection 2.2). Neglecting contributions from radii smaller 
than a cutoff radius R the reduced radial integrals can be approximated by 



cos{5ij)T\l{X)+sin{5i^)g\l{X) 



(47) 



where (A) and Q\^. (A) denote the real and the imaginary part of the function 



^!f(A)=g"+^ / drr' Fi^{r]f;kr) + zGi^{r]f,kr) (2gr) (48) 



,A+1 



R 
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j dtt^ [Fi^{rii/x;x^t)+iGi^{rii/x;x^t)\ W_^^^i^^{2^t) 



Here, as in the following, the quantum numbers Jj, ji, Jf, jj and have 
been suppressed if no confusion arises. The nuclear phase shifts encode 
the effects of the final-state interaction. For 5;^ = one obtains the results 
without the nuclear interaction between the nucleon and the core. Note that 
the function 7i/. (A) depends only on the parameters 7, rji and x. 

If there is only one fixed pair (jj, k) in the initial state and similar {jf, //) in 
the final state the expression for the reduced transition probability reduces to 



dE 



{EX, Jisjc kJfsjc) = 



(49) 



2 2// 2Jf + 1 
nh^ 2Ji + 1 



Ji 1 



2A+3 



with the dimensionless shape function of the transition strength 



(50) 



that completely contains the dependence on the momentum in the contin- 
uum. Similarly, one can define the characteristic shape functions for the photo 
absorption 

5gi^^(abs, A,.) = (1 + x^f-' sil^^{Xj^ (51) 

and for the capture cross section 

S^:f^, (capt, Aje) = ^ ^ J S2^£^ ( Aje) . (52) 

The expression (49) directly shows the dependence of the transition strength 
on the characteristic parameters Cj.j^i^ and q of the ground state. For small 
separation energies Sb of the nucleon the strength becomes very large since q is 
a small number. Note that the ANC in general depends on q, too. In a model 
for neutrons in a square- well potential one finds a dependence C-j.^.i. oc for 
k — O and Cj^j oc g'* for 1^ > 0, respectively, for small q (see Appendix A) . 

When the nucleon 6 is a neutron the reduced radial integrals can be calculated 
analytically (see Appendix C). Then the functions (47) only depend on the 
dimensionless variables 7 = qR and k — kR — x^ and the phase shift in 
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Table 2 

Characteristic shape functions S^^ (A) for 7 = and no final-state interaction. 



A — U 




A = 1 




A = 2 




50(0) = 

5^(0) = 


X 


5?(1) = 

50^1) = 
<52'(1) = 
5f (1) = 


x(3 + x'^y 


50(2) = 

Sl{2) = 


x(15 + 10^2 + 3x4)^ 


(1+X2)2 
X^ 


(l + a;2)4 
4x^ 


(l + a;2)6 
4x^(5 + x2)2 


(1+X2)2 


(1 + .t2)4 

^■■^(S + 3a;2)2 


(l + x2)6 






(1+X2)4 


5^(2) = 


64x5 






(1+X2)4 


(1+X2)6 



the final state. It is found that the reduced radial integrals have the general 
form 



X 



7 e xp(-7) 

nlf'^ (A) cos(«: + 5, J + TeS")'^ (A) sin(/. + 5,,; 



with polynomials/rational functions TZ[f^''^ {X). Explicit expressions are given 
in Appendix C. In general they are complicated functions of 7 and k. 

From the discussion in subsection 2.4 a x~'^ dependence of the reduced radial 
integrals X^'/(l) for dipole transitions is expected for large x. However, the 
integrals (53) for A = 1 show a different high-a; behaviour and the convergence 
of the sum rules (see section 4) is not guaranteed. This is a consequence of 
neglecting the interior contribution to the radial integral that becomes relevant 
at high relative energies since the nucleon penetrates into the core. Thus, 
the reduced radial integrals (53) are only a good approximation for not too 
high relative energies in the continuum. In comparison, equation (29) is valid 
irrespective of the halo nature of the system and for all energies but only 
for A = 1 and the square-well case; it contains all contributions to the radial 
integral from zero to infinity. Equation (53) is a good approximation for halo 
nuclei at small relative energies, since the interior contributions are small in 
this case (cf. Fig. 4). It can also be applied easily to cases where the potential 
is different in the bound and scattering states. It is worthwhile to study some 
limiting cases. 



3.1 Shape functions in n+core systems without final- state interaction 



Without nuclear interaction between the neutron and the core in the final state 
the phase shift Sij, is zero and there is no contribution from the integral with 
the irregular wave function in (47). The contribution for radii r < it! in the 
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Fig. 6. Dependence of the generic shape functions S^^ (X) with cutoff radius R = 
on = E/Sn for various transitions li If of multipolarity E1/E2 (top/bottom). 
The functions are normahzed to 1 at = 1. 

radial integral is small and it is possible to take the hmit R ^ keeping the 
neutron separation energy Sn or equivalently the inverse bound-state decay 
length q constant. In this limit both 7 and n approach zero but the ratio 
X = k/j = \J E / Sn, i.e. the relevant variable for the shape of the strength 
distribution, is independent of R. In this case the reduced radial integrals 
(53) assume a particular simple form. They only depend on this dimensionless 
variable x. Prom the reduced radial integrals the shape functions (50) of the 
reduced transition probability are easily derived. The functions (X) in the 
limit R with 61^ = are given in table 2. Especially the expression iSq(1) 
is well known, see, e.g., [13,45]. It can be found in a different notation also in 
Refs. [2,37]. 

In Fig. 6 the generic form of the shape functions Sj^ (A) is shown as a function 

of x'^ = E / Sn for El and E2 transitions from a bound state with orbital angu- 
lar momentum /j to a scattering state with orbital angular momentum If. The 
dependence of the transition shape on the centrifugal barrier is clearly seen. 
The peak at small energies is more pronounced for low orbital angular mo- 
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Fig. 7. Shape functions (X) as a function of = Ehc/ Sn for values of the 
parameter 7 between 0.0 and 1.0 in steps of 0.2. 

menta //. Only for s and p waves in the continuum a large transition strength 
is found close to the threshold. At low x we have the typical dependence 

5lf(A)ocx^'/+^ (54) 



The shape function for photo absorption (abs, A) has the same x depen- 
dence as (A) for small final state momenta. The shape function for radiative 
capture 5;*(capt, A), on the other hand, shows a x'^''f'^ dependence for small 
momenta in the continuum due to the additional momentum dependent factor 
from the theorem of detailed balance (40), see equation (51). 

In case of a finite cutoff radius R in the radial integral (18) the shape function 
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(A) depends on both 7 = qR and k — kR. Since q is fixed for a given 
nucleus is it reasonable again to use 7 and x = n/'j = k/q as independent 
parameters. For halo nuclei with small nucleon separation energy q will be a 
small quantity. As long as R does not become too large, e.g. for heavy nuclei, 
7 will also be small. The variation of the shape functions with 7 gives an 
estimate of the contribution to the radial integral from the nuclear interior. In 
Fig. 7 the change of the shape functions with 7 is shown for the transitions of 
Fig. 6. There is a clear systematic trend. The shape functions are less sensitive 
to a change in 7 for larger final state orbital angular momentum If and higher 
multipolarity A because of the suppression of the integrand at small r from 
the spherical Bessel functions jif{kr) and from the transition operator, 
respectively. On the other hand, a larger orbital angular momentum Zj in the 
bound state increases the sensitivity since the wave function introduces a r~^' 
dependence at small radii. This explains the strong 7-dependence of the shape 
function 52(1). 

3.2 Shape functions in n+core systems with final- state interaction 
For neutron scattering the finite-range expansion (41) reduces to 



since r/ = 0. Taking only the contributions with ai and ri into account the phase 
shift Si crosses the value 7r/2 at an energy Eq = h'^/{fiairi). This behaviour 
produces a resonance in the corresponding partial wave if the phase shift is 
increasing that, however, does not have to be physical. In general, the actual 
energy-dependence of the phase shift close to a resonance at energy E^i is given 
by the Brcit-Wigner form tandi = V / {2{Er — E)) with positive width F where 
the parameters E^i and F are not related to a/ and in the effective-range 
expansion. Particular values for these parameters that correctly reproduce the 
phase shift at low momenta will not necessarily reproduce the position and 
the width of an actual resonance with orbital angular momentum I. However, 
the phase shift can be calculated with the help of the relation 



for general cases assuming an energy dependence of the function hi. It replaces 
the phase shift 5i in order to take the FSI in the scattering wave function into 
account. The dimensionless function hi is related by 



0+1 




(55) 




(56) 




21+1 



(57) 
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to the function Di{E) introduced in Ref. [46] for a low-energy expansion of 
cross sections for radiative capture reactions. 

The function hi is of the order ol ^ — qR unless the logarithmic derivative 
L of the scattering wave function at radius it! is close to —I, see eq. (A. 31). 
It is useful in expansions of the shape functions (A) if the limit — > is 
considered. On the other hand, the scaled function 



(usually of order one) is the appropriate quantity if the limit 7 ^ is studied. 
In this case we have 

ai = lim {ciRf^^ (59) 



with constant radius R. 

Typical values of hi and q can be estimated from examples employing a single 
particle model where a nuclear potential of Woods-Saxon form with typical 
parameters is assumed. In Figure 8 the function hi for orbital angular momenta 
/ = 0, 1, 2 is shown as a function of x = k/q for the scattering of a neutron on 
^''Bc assuming different depths Vq of the potential. In most cases the function hi 
is quite constant as a function of x except when there is a resonance (vertical 
lines in Fig. 8) in the continuum for a certain fixed depth Vq. For p and d 
waves hi is usually in the intervall [—0.5,0.5]; for s waves 60 covers a larger 
range. Assuming a zero-range potential for the neutron-core interaction the 
scattering wave function is given by ijj^'^\r) = exp{ik-r) — exp{ikr) / [{q + ik)r] 
with the bound state parameter q determined by the binding energy. The s 
wave scattering length is just gq = 1/q corresponding to 60 = 1 and cq = I/7. 
Since the ground state is close to the threshold for a halo nucleus cq becomes 
unnaturally large and there will be a large effect on the continuum in the 
same partial wave. However, the partial waves in the final state for an electric 
excitation have a different orbital angular momentum and parity and usually 
there will be no resonance in the energy range of interest. In general it is 
reasonable to assume < 0.5 with a weak momentum dependence. 

The phase shift in the general expression for the reduced radial integral (53) 
can be replaced by the function 6; using the explicit relation (56). For some of 
the transitions it is still possible to take the hmit — > in the reduced radial 
integrals for hi ^ 0. The corresponding shape functions are given in table 3. 
For the other transitions the reduced radial integral and the shape function 
diverges. It is obvious that the x^''^'^^ dependence at small x is obtained again 
but the modulus is modified in the limit x — > as compared to the generic 
shape functions in table 2 because of the FSI. 
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Fig. 8. Function bi in partial waves with orbital angular momentum / = 0, 1,2 as 
a function of x = k/q for the scattering of a neutron on ^''Be in a single-particle 
model with a Woods-Saxon potential of radius Ro = 2.78 fm, diffuseness parameter 
a = 0.65 fm and various potential depths Vq. (Note the nonlinear scale on the y 
axis.) 

In Fig. 9 the dependence of the shape function in the limit i? — > on 6/^ is 
shown for the cases of table 3. (Transitions with A = are not relevant.) The 
function varies in the interval [—0.5,0.5]. This corresponds to reasonable 
values that are far from a resonance in the particular continuum channel. 
Depending on the magnitude of bij, the shape functions show a pronounced 
variation around the generic form that reflects the effect of the potential in the 
final state. Small changes in the reduced scattering length lead to a smaller 
effect for higher orbital angular momenta If due to the bf -^'^^ dependence in 
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Table 3 

Characteristic shape functions 5^^:' (A) for 7 = and bi^ / 0. 



A = 




A = 1 




<5o°(0) = 


X (1 - feo)' 


<S?(1) = 
5o'(l) = 


X (3 - 26o + ^2)2 


(1+X2)2 l + (x6oP 


(l + x2)4 l + (a;6o)2 
x3 [2-6?(l + 3x2)]2 






(1 + ^2)4 l + (a;6i)6 


A = 2 



5°(2) 



(1 + X2)6 1 + (x5o)2 

1. . ^ x^ 4[5 + a;2-fe3(i + 5x2)]2 

1^ > (1+^2)6 1 + (x6i)6 

5 [8-6|(3 + 10a;2 + 15x4)]2 



°^ ^ (1+X2)6 1 + (X62) 



10 




Fig. 9. Shape functions {X) for 7 = as a function of x^ = E/S^ for various 
values of the function bij, in steps of 0.1. The thick dashed, solid, and dot-dashed 
lines correspond to 6;^ = +0.5, 0.0, —0.5, respectively. 
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Table 4 

Expansion of the characteristic shape functions S^^ (X) in the parameter x for finite 
7 and 6/^ / 0. 

A = 

50(0) = exp(-27) (1 + 7- 6o)'x 

5,^(0) = g^^^^ [7(3 + 37 + 7^) - Sbl] ' x' 

A = 1 

5iO(l) = exp(-27) [3 + 37 + 7^ - bo{2 + ^)]\ 

ci/iA exp(-27) . , Q 2 , 3 o/,3\2 3 

■SoW = g (6 + 67 + 37 +7 -361) X 

'^2(1) = '-^^^ [^(15 + + 67^ + 7^) - 36?(3 + 7)]^^ 
<g?(l) = [7(30 + 3O7 + 457^ + 57^ + 7^) - 456i] ' x' 

A = 2 

520(2) = exp(-27) [l5 + I57 + 67^ + 7^ - 60(8 + 57 + 7^)]% 
51(2) = !^pLM [30 + 3O7 + 157' + 573 + 7' - 36?(2 + 7)] ' x^ 
S^{2) = ^^^4=^ (120 + I2O7 + 6O72 + 2O73 + 57^ + 7^ - 456i)' 



the analytical expressions. For positive values of bij_. one finds an increase of 

{X) whereas a negative bi^. leads to a decrease in absolute value. There is 
also a change in the shape observed and a shift of the maximum. This shift 
becomes noticable only for large values of 6;. The position of the maximum 
moves to larger x"^ with decreasing hi for p and d waves in the final state, 
whereas the trend is opposite for s-wave final states. There is also a clear 
hierachy observed. The maximum of the shape function appears at higher 
with larger orbital angular momentum If in the continuum. 



The analytic expressions of the shape functions in table 3 were obtained by 
extending the asymptotic form of the bound and scattering wave functions to 
zero radius in the radial integral. The integrand with the regular scattering 
wave function shows a r'^"''"'"'^"'"^ dependence; in contrast to that one finds a 
^x-h-if behaviour at small r for the integrand with the irregular scattering 
wave function. One might expect that the irregular contribution is overesti- 
mated in this limit since it diverges for r — > 0. Expanding the general shape 
functions in powers of x with constant 7 — qR one finds the results as given 
in table 4. The functions again show the typical x^'-f"*"^ dependence at low 
X = k/q but the slope depends less strongly on bi^. for finite 7 than for the 
functions given in table 3. The shape functions of the transitions d ^ p and 
p — > d (p — > p) for A = 1 (A = 0) diverge in the hmit 7 — > or i? — > 0. 
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Fig. 10. Shape functions (X) for 7 = 0.5 as a function of = E/Sn for various 
values of bij, in steps of 0.1. The thick dashed, sohd, and dot-dashed hnes correspond 
to bij: = +0.5, 0.0, and —0.5, respectively. 



Fig. 10 shows the dependence of the shape functions on bi^ for constant 
7 = 0.5, i.e. a constant cutoff radius R. This value of 7 corresponds to 
R = 3.36 fm in case of neutron scattering on ^'^Be with a bound state pa- 
rameter q = 0.1487 fm~^ for a neutron separation energy of Sn = 0.504 MeV. 
The shape functions are very similar to the case with i? — > 0, cf. Fig. 9, except 
for two cases where (A) diverges in the hmit — > for finite bi^. In general 
the shape functions for finite R are slightly smaller than in the case R — 0. 
Final state effects for a finite value of 7 are less pronounced than in the case 
with 7 = but the differences is small. Therefore the dependence of the shape 
functions in table 3 on bi^ gives a reasonable impression about the importance 
of the final state interaction. For El transitions p— >dord— >pa strong 
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2 

X 



Fig. 11. Shape functions (X) for 7 = 1.0 as a function of = E/Sn for various 
values of bij, in steps of 0.2. The thick dashed, solid, and dot-dashed lines correspond 
to bij: = +1.0, 0.0, and —1.0, respectively. The scale of the y axis is chosen to be 
the same as in Fig. 10 for a better comparison. 



dependence of the shape function on bij_. is found. The shape and magnitude 
of the reduced transition probability depends sensitively on the choice of 6/^ 
and R. 

The sensitivity of the shape functions to final-state effects on the neutron 

separation energy can be estimated from an expansion of 5/ (A) in terms of 
7 by replacing the function hi with the quantity q as defined in equation 
(58). This approach was introduced in [25]. Analytical expressions for dipole 
transitions are given by 



32 



a;(3 + x' 
Ax^ 



2\2 r 



1 - 



4co 
3 + x' 



r7 + 



x3(5 + 3a;2)2 



:i + a;2)4 
Ax^ 



[l + x 



2\A 



1 - 



5 + 3x2 ^ +■■ 



60 



(60) 
(61) 
(62) 
(63) 



These expansions are consistent with the results in table 4 in the limit 7 — > 
(q = 6;/7, see eq. (58)). There are corrections to the El shape functions in 
table 3 linear in 7 for p — > s transitions and proportional to 7^ for s ^ p 
transitions for finite values of c/^. For d — > p and p — > d transitions one 
finds corrections to the forms in table 3 of order 7^ and 7^, respectively. They 
contain a contribution depending on the interaction (q^,) and a correction from 
the finite size of the system. Transitions / — > / — 1 are affected more strongly 
by the interaction than transitions I ^ I 

Another approach to estimate the sensitivity to the FSI is a comparison of 
the shape functions by scaling 7 and 6/ appropriately. In Fig. 10 the shape 
functions were shown for 7 = 0.5 and constant hi^ G [—0.5, 0.5]. Multiplying 7 
by a factor of 2 corresponds to an increase of S'„ by a factor of 4. In order to 
cover the same range of scattering lengths ai as in Fig. 10 also the interval for 
the function hi has to be increased to [—1.0, 1.0] because the ratio hi^/q has to 
be kept constant, cf. Eq. (45). The corresponding shape functions for 7 = 1.0 
are depicted in Fig. 11. It is obvious that a higher neutron separation energy 
leads to a much stronger dependence of the shape functions on the strength of 
the final-state interaction. This behaviour is easily understood from the radial 
dependence of the bound state wave function. For larger separation energy 
the slope of the asymptotic wave function is much steeper since q increases. 
The radial integral is more sensitive to a small shift of the continuum wave 
function from a finite nuclear phase shift because the main contribution to 
the integrals arises from a smaller intervall in the radius on the surface of 
the nucleus. The effect is less pronounced for s waves in the final state but 
increases dramatically for higher partial waves in the continuum. 



3.3 Shape functions in p+core systems 



The Coulomb interaction in p+core systems leads to a systematic modification 

of the characteristic shape functions as compared to n+core systems. Since 
analytical results are not available in general one has to resort to a numerical 
integration of the radial integral (48) if the dependence on arbitrary values of 
X is studied. In Ref. [47] an analytic result for the radial integral was obtained 



33 




Fig. 12. Scaled shape functions e^*5^:'(A) as a function of = Ebc/Sp for 7 = 
and values of the parameter rji between 0.0 and 1.0 in steps of 0.2 without nuclear 
final-state interaction. 

in the special case of scattering on a solid sphere with given radius. For x ^ 1 
it is possible to obtain analytical results in an approach with an expansion for 
small energies as presented in Refs. [41,48]. For small x one finds a suppression 
of the shape function approximately proportional to exp(— 27r?7j) with the 
Sommerfeld parameter Tjf of the scattering state. This scaling is characteristic 
for a case with a Coulomb barrier. However, at larger energies the absolute 
value of the shape functions is determined by the Sommerfeld parameter rji of 
the bound state since it defines the range of radii with the largest contribution 
to the radial integral. In figure 12 the variation of the shape functions with r/j 
for 7 = and for various transitions is shown. The functions are scaled with 
exp(4?7j). Effects from the nuclear interaction in the final state are not taken 
into account. For rji — the results for neutrons are recovered. With increasing 
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Fig. 13. Ratios R/.{X,r]i) function of = Ehc/ Sp for 7 = and values of 
the parameter rji between 0.0 and 1.0 in steps of 0.2 without nuclear final-state 
interaction. 



r]i the maximum of the shape function shifts to higher x^. i.e. to higher relative 
energies. At the same time the width of the strength distribution increases and 
the absolute value of {X) reduces considerably (note the scaling). 

At low x"^ the suppression of the shape functions with finite rji as compared 
to the case rji = is quantitatively given by the factor Cf^{rif) as defined 
in equation (43). For x ^ 00 corresponding to r^/ = rji/x — > we have 
hm^^_>oCi^(?7/) — 1. This suggests to scale the shape functions for the pro- 
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ton+core case with the penetrabihty factor Ci '^{r}f). In figure 13 the ratio 



2 Si!i\vi) 



(64) 



is depicted for various transitions and parameters rji. The ratio only weakly 
depends on x'^. Therefore, the main difference between the proton+core and 
neutron+core cases is well described by the penetrability factor C^~^(?7/). 

There is a major difference between n+core and p+core nuclei. Due to the 
increasing Coulomb barrier the observation of a large transition strength at 
low relative energies will be lost for heavier p+core nuclei and the difference 
between transitions of different orbital angular momenta in the initial and 
final states will be become smaller. In contrast, even for heavy n+core nuclei 
the strong transition strength at small energies will persist and the shape will 
still be characteristic of the particular transition. In principle, variations of the 
generic shape functions as shown in figure 12 with finite parameters 7 = qR 
and bi can be studied along the lines as for the n+core case. Very similar 
trends are observed and we omit a detailed discussion. 



4 Total transition strength and sum rules 

The total strength for an EX transition from a bound state with quantum 
numbers Ji,s,jc,li to all possible final states is related by the non energy- 
weighted sum rule 



to the expectation value {r'^'^)i^ of the initial bound state with orbital angular 
momentum /j, see, e.g., [49]. Note that for the integration over the energy 
E not only continuum states but also bound states have to be taken into 
account if they can be reached by an EX transition from the initial bound 
state. In the special case of dipole transitions the root-mean-square radius of 
the bound state wave function determines the total El strength. From the 
scaling laws of {r'^)i^ as discussed in subsection 2.2 we expect a divergence of 
the total transition strengths B{El,li) oc 7"^ and B{El,li) oc 7"^ for li = 
and 1, respectively, in the hmit 7 — > 0. In contrast, the total El strength 



B{EX, JiSjJi) 



j dE -^{EX, Jisjck kJfsjJf) 

Jfh 



(65) 
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should remain finite for larger orbital angular momenta of the bound state 
wave function. 

The total reduced transition probability to the continuum is obtained in our 
approach with asymptotic wave functions from 



Bcont{EX, JiSjJi JfSjJf) 



(66) 



2^/ + l 
2Ji + l 



^2A+l 



with the dimensionless integral 



(67) 



Evidently, the continuum contribution is only a fraction of the total strength if 
the FSI allows the occurrence of bound states in the corresponding channels. 
Since wc use only the exterior contribution to the radial integral it is not 
guaranteed that the high energy behaviour of the reduced radial integrals 
(53) is correct. Indeed, generally we do not find the x~'^ dependence for El 
transitions as expected from the general considerations in section 2.5. For large 
X there is a sizable contribution from the interior integral that is necessary 
to generate the correct x dependence at high energies. However, the main 
contribution to the integral (67) arises at small x close to the threshold and 
as long as the integral converges it will give a reasonable approximation of the 
total transition strength. 

If one assumes that there is no nucleon-core interaction in the final states (i.e. 
plane waves) the full strength lies in the continuum. In this case the func- 
tions Ti_^ (A) can be calculated analytically for neutron+core systems in our 
approach. These functions are given in table 5 where we suppressed unim- 
portant quantum numbers for simplification. The integrals ^^(0) and ^^(1) 
diverge in the limit 7^0 but remain finite for 7 > 0. In the limit 7 —> 00, i.e. 
if the bound state does not develop a halo, the functions become very small. 

The total transition strength (66) is proportional to |C;. |^ 7^'^(A)/g^'^+^ The 
scahng of the ANC Cl^ in the limit 7 — > can be obtained by considering the 
square- well model (see appendix A). We find 



q if = 

(68) 

572''-! iili>l 
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Table 5 

Functions 7^,^ (A) for finite 7 and 6;^ = 0. 
A = 



ro°(0)=exp(-27)/2 

rii(0) = (2 + 7)exp(-27)/(27) 



A = 1 




= (5 + 67 + 27^) exp(-27) /4 




= (l + 27 + 72)cxp(-27)/4 




= (36 + 377 + 147^ + 27'^) cxp(-27)/(47) 


r2(i) 


= (5 + 67 + 27^) exp(-27)/4 


A = 2 




= (63 + 9O7 + 547^ + 167^ + 27'') exp(-27)/4 




= (7 + 147 + 147^ + 87^ + 27^) exp(-27) /4 


^o'(2) 


= (3 + 67 + 67^ + 47^ + 27^) exp(-27) /4 



for 7—^0. Taking this dependence into consideration we find the scahng of 
B{El,li) consistent with the expectation from the non energy- weighted sum 
rule. 

Now, let us consider the case with non-vanishing nuclear interaction in the 
continuum states. The FSI is parametrized in our approach by the function 
bi^ that in general depends on the energy E. If we assume that bi^ is constant 

for all energies the functions if^iX) can be calculated analytically again. One 
obtains, e.g. for the El p ^ s transition 



" 5+67+27" _ 2/. [i+2feo+7(i+feo)]" (-^]] for 6n > 



5+67+27-^ 



for bo<0 



For other transitions the expressions become quite complicated and we abstain 
from giving them here. Since the function 6/^ cannot be assumed constant over 
the whole range of continuum energies one cannot expect that the calculated 
{X) corresponds to a reasonable result in general. However, it is sufficient 
to assume that the function bi^ does not vary too much over the peak of the 

transition strength that is described by the shape function s\^. (A) . For large 
I fool the excitation spectrum is distorted significantly in the relevant low-energy 
region. 

Equation (69) is quite instructive. It gives a good idea how the B{E1) strength 
is modified when the depth Vq of the nuclear potential (with reasonable ge- 
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Fig. 14. T^{1) as a function of 1/bo for values of the parameter 7 between 0.0 and 
1.0 in steps of 0.2. 

ometry) changes. For 60 > there is a reduction of the total strength in the 
continuum whereas in the case 60 < no reduction is found. This behaviour 
is related to the occurrence of bound states in the s-wave continuum. For 
Vo = the nuclear phase shift Sq and the function bg is zero. Increasing the 
potential depth the phase shift at low energies becomes positive, correspond- 
ing to a negative scattering length. Increasing the depth Vq further will lead 
to a more negative 60 and a resonance will appear at low energies. Then there 
comes the point where the potential will be able to support a bound state and 
60 suddenly jumps to a large positive value. Increasing the depth Vq further 
will reduce the scattering length again. Finally, 60 becomes zero and the cycle 
starts again. Considering this relation between strength of the potential and 
the scattering length it is reasonable to plot as a function of I/60 as 

depicted in figure 14. As long as the potential is not able to bind a state the 
continuum B{E1) strength does not change. As soon as a resonance in the 
continuum becomes bound a sudden drop of the transition strength is observed 
that recovers to its initial value when approaches zero from positive values. 
An analogous effect was found for the s ^ p transition for the El excitation 
of ^^Be in Ref. [25]. Similar observations are expected for other transitions. 

This effect of the FSl has consequences for the extraction of the ANC or 
the spectroscopic factor for an assumed ground-state single-particle config- 
uration from experimentally obtained total reduced transition probabilities. 
When the value of B{E1, li) from the experiment is compared to the theoreti- 
cal result assuming a single-particle model with plane waves in the continuum 
the extracted ANC or spectroscopic factor will be underestimated if the ac- 
tual nuclear potential supports also bound states. In the latter part of 
the total transition strength goes to these bound states, even if they cannot 
be reached by a transition of the halo nucleon because they are occupied by 
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nucleons of the core. 

Using the relation 

j:iEf-E.)\{mi)f^l{i\H[H,r,]]\i) (70) 
the energy-weighted (or Thomas-Reiche-Kuhn) sum rule 

oo 

S{E1, JiSjJ,) = ^Y.JdE{E + S,) —{El, J.sjM ^ kJfSjJf) (71) 
■Jfh 

- 47r r eff 

for dipolc transitions is obtained. In constrast to the non energy-weighted 
sum rule (65) it gives a result independent of the orbital angular momentum 
li of the bound state. In principle, the sum contains all bound and unbound 
states that can be reached from the initial state by an El transition. Since the 
photon energy is E^ — Ef — E^ — E + S the sum rule (71) can be expressed 
in terms of the total El photo absorption cross section (38) as 



j dE^ (TEi{a + 'y ^ b + c) 



(72) 



2 27r2/le2 
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N,Z, 



NcZ, 



for the disintegration into core-|-nucleon. It has the form of a cluster sum rule 
since we only consider the excitation of the halo nucleon but not of the core 
[50,51]. In our approach we can write for the continuum contribution 



'S'cont 

{El, JiSjck) 



(73) 



2Ji + l 



with the dimensionless integral 
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If there arc bound states in the final channel we expect S'cont (-£■!, Jisjch) < 
S{E1, JiSjck). Here again, the remarks on the high energy behaviour of the 
integrand apply as for the functions (67). For the most important transition 
s — > p without interaction in the continuum states we obtain the finite result 

Wo^(l) = e--^ (75) 

where we suppressed irrelevant quantum numbers. Considering the scaling 
law (68) of the ANC for li — we find the correct result (71) for the energy- 
weighted sum rule in the limit 7 — >^ 0. For the ground state with k = 1 there 
are two contributions from // = and // = 2 in the final state. In this case 
the integral (1) diverges with 7"^ but this dependence is compensated by 
the scahng of the ANC. However, since the radial integrals in the external 
approximation do not vanish fast enough with increasing x, the correct value 
for the energy- weighted sum rule is not recovered in the limit 7 — > 0. 



5 Examples for transition strengths of nucleon+core nuclei 

Let us now look at some specific examples of exotic nuclei with nucleon+core 
structure. The dependence of the reduced transition probability and of the 
cross sections for electromagnetic transitions on the strength of the fragment- 
fragment interaction in the final state can be studied in the simple but realistic 
single-particle model of subsection 2.3. For the nuclear potential a Woods- 
Saxon shape with the representative parameters a = 0.65 fm and Rq = TqA^^^ 
where Tq = 1.25 fm is assumed. The depth of the potential in the bound state 
is adjusted to reproduce the experimental neutron and proton separation en- 
ergies, respectively. In order to simulate a varying interaction strength in the 
final states the transition strength is studied for the depth of the Woods-Saxon 
potential in the continuum in the range from to 80 MeV. In Tables 6 and 7 
the separation energy Sb, potential radius Rq, potential depth Vq, and orbital 
angular momentum /j for the ground state are given for the neutron+core 
nuclei and proton-|-core nuclei, respectively. For simplicity, no spin-orbit po- 
tential is considered in the calculations. Also the characteristic parameter 7 
is given in Tables 6 and 7. Only nuclei with 7 < 1 can be considered to be 
halo nuclei, however, the boundary to nuclei without an extended nucleon dis- 
tribution is smooth. For p-l-core nuclei the parameter rji is given in addition. 
A larger value indicates a stronger dominance of Coulomb distortions of the 
dB{El)/dE strength. In the case of neutron-|-core nuclei £'1 s — > p and d ^ p 
transition will be considered. In the case of proton+core nuclei we will study 
especially El p s excitations. For many observables, a square- well model 
would also be a good approximation, especially for low energy processes where 
the "shape independence" is valid to a large extent. 
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Tabic 6 

Neutron separation energy Sn, radius Rq, depth Vq of the Woods-Saxon potential 
and parameter 7 = qRq for the ground state with orbital angular momentum Zj of 
the neutron+core nuclei in the single-particle model. The diffuseness parameter is 
a = 0.65 frn^ 
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0.413 
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5.1 Neutron+core nuclei 

The 1^ ground state of ^^Bc is considered as a neutron in a 2si/2 state cou- 
pled to the 0+ ground state of ^°Be. The experimental neutron separation 
energy of Sn = 0.504 MeV is quite small so that ^^Be is a typical example 
of a halo nucleus. It was studied by Coulomb dissociation in the experiments 
[45,52,53,54]. In Fig. 15 the reduced transition probability dB{El) / dE for 
a transition into p waves is shown for varying depths Vq in the continuum 
between and 80 MeV in steps of 5 MeV. The El strength shows a large 
peak at relative energies well below 1 MeV typical for a halo nucleus. The 
shape changes smoothly with increasing potential depth except for Vq close to 
30 MeV where a p-wave resonance in the continuum becomes bound. For a 
potential depth similar to the one required for the correct separation energy 
in the ground state (cf. Tab. 6) the shape of the El transition strength is very 
similar to the result for a plane wave, i.e. Vq = MeV. Experimental data for 
the dB{El)/dE distribution in ^^Be were described in [25] in our approach 

by fitting the asymptotic normalization coefficient and the constants c^^^ and 
1/2 

Ci in the effective-range expansion for the phase shifts in the p waves with 
total angular momentum j = 3/2 and j = 1/2. An unnaturally large value 
for Ci^^ was found that is related to the low separation energy of the first 
bound excited state in ^^Be. The shape of the dipole strength function in the 
effective-range approach was found to agree very well to a calculation in the 
Woods-Saxon model. A corresponding figure and more details can be found 
in Ref. [25]. 

Another example of a neutron-|-core nucleus with the neutron in a s wave 
ground state is ^^C but with a larger neutron separation energy of Sn = 
1.218 MeV [55]. Comparing the dependence of the dipole transition strength 
dB{El)/dE on the potential depth Vq as shown in Fig. 16 with the case of 
^^Be one observes a stronger variation of the shape and absolute value. This 
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dB(El)/dE 
[e^fmVMeV] 




Fig. 15. Reduced transition probability dB{El) / dE for the breakup of ^^Be into a 
neutron and ^°Be as a function of the cm. energy E for various depths Vq of the 
potential in the contiuum. The quantum numbers for the single particle bound state 
and the neutron separation energy in MeV are given in parenthesis. 

clearly shows the increased sensitivity of the transition strength on the final- 
state interaction when the neutron separation energy increases as expected 
from the analytical model in subsection 3.2. Again, the reduced transition 
probability changes quite abruptly when a p wave resonance in the continuum 
becomes a bound state. 

An even more drastic change of the transition strength is observed in the case 
of ^•^O when the neutron is assumed to be in a s-wave ground state, see Fig. 17. 
This nucleus with a neutron separation energy of Sn — 2.74 MeV cannot be 
considered as a real halo nucleus. Yet, it can be said that there is substantial 
low-energy El strength. 

Comparing the three cases ^^Be, ^^C, and ^^O one observes a shift of the 
peak in the transition strength to higher relative energies in the continuum. 

The position of the peak scales with the separation energy as expected from 
the dependence of the shape functions on the ratio = Sn/E in Section 3. 
Additionally, one finds a reduction of the overall strength as suggested from 
Eqs. (49) and (66) with an increasing value of the ground state parameter q. 

Since the ground state angular momentum of ^^O is not known uniquely from 
experiments [26,56,57,58,59,60,61] it is also possible that the neutron in the 
ground state occupies a 0^5/2 single-particle state. The strength function for 
an El transition to p-wave final states under this assumption is shown in 
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Fig. 17. Same as Fig. 15 but for the breakup of O in a s-wave ground state into 
a neutron and ^^O. 

Fig. 18. The absolute magnitude is about a factor 4 smaller than in the case 
of a s-wave neutron in the ground state and the maximum is shifted to larger 
relative energy. This is explained by the additional centrifugal barrier in the 
d wave that reduces the probability of finding the neutron at large distances 
from the core in the classically forbidden region. 
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O (ld5/2, 2.74 MeV) 



dB(El)/dE 
[e^fmVMeV] 




Fig. 18. Same as Fig. 15 but for the breakup of O in a d-wave ground state into 
a neutron and ^^O. 



O (ld5/2, 4.14 MeV) 



dB(El)/dE 
[e^fmVMeV] 




Fig. 19. Same as Fig. 15 but for the breakup of "O into a neutron and °0. 

A further example for a neutron+core nucleus with a o?-wave ground-state 
configuration is ^^O, a nucleus close to the valley of stability. Here the sepa- 
ration energy of S'„ = 4.14 MeV is even larger than in the case of ^'^O and the 
shape and magnitude of the reduced El transition probability to p waves in 
the final state vary drastically with an increasing strength of the final-state 
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Fig. 20. Total reduced transition probability B{E1) for the breakup into a neu- 
tron+core with cm. energies E between and 10 MeV as a function of the potential 
depth Vo for the continuum states. The filled circle gives the result for the same 
potential depth as for the single particle bound state. The quantum numbers for 
the single particle bound state and the neutron separation energy for each system 
are given in parenthesis. 

interaction. 



In Fig. 20 the total El transition strength integrated from to 10 MeV is 
shown for the nuclei discussed above as a function of the potential depth 
Vo in the continuum. There is a distinctive difference between nuclei with 
a s-wave and li-wave neutron in the ground state. In the former case the 
integrated strength stays almost constant when Vq is increased starting at 
MeV, i.e. the plane-wave result. Then there is a sudden drop of the total 
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Table 7 

Proton separation energy Sp, radius Rq, depth Vq of the Woods-Saxon potential, 
parameters 7 = qRo and rji for the ground state with orbital angular momentum Zj 
of the proton+core nuclei in the single-particle model. The diffuseness parameter is 
a = 0.65 fm. 











1 2 AT 


Sp 


[MeV] 


0.1375 


1.296 


0.601 


Rq 


[fm] 


2.50 


2.60 


2.86 


Vo 


[MeV] 


43.183 


44.4456 


36.474 


li 




1 


1 


1 


7 




0.190 


0.613 


0.466 


Vi 




1.595 


0.655 


1.171 



El strength when the p-wme resonance from the continuum becomes a bound 
state. However, the total strength recovers to its value for = when the 
potential depth is increased further until the next higher p-wave resonance 
becomes bound. The dependence of the total transition probability on the 
strength of the FSI is in qualitative agreement with the expectations from the 
analytical neutron-|-core model as shown in figure 14. The rise of the total 
El strength after the continuum to bound state crossing is faster when the 
neutron separation energy is smaller. For halo nuclei there is only a small 
variation of the total reduced transition probability if one is not too close to 
a potential strength where the crossing appears. When the neutron is in the 
d-wave ground state there is already a considerable dependence of the total El 
strength on the depth of the potential close to the plane- wave limit. Beyond 
the sudden drop of the strength at the continuum-to-bound-state crossing the 
absolute value increases again but it does not reach the plane-wave limit again. 
Comparing the total El strength calculated for a continuum potential that is 
identical to the bound state potential (filled circle in Fig. 20) one immediately 
finds that the integrated El strength is almost the same as in the plane-wave 
case for a s-wave neutron in the ground state but it is significantly smaller for 
a d-wave neutron. 



5.2 Proton+core nuclei 

Systems with proton-|-core structure show similar features as compared to 

neutron+core systems, however, with modifications due to the appearence of 
the Coulomb barrier. Since analytical results for the shape functions are not 
available one has to resort to the numerical calculations, e.g. in the simple 
single-particle model. 
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Fig. 21. Reduced transition probability dB{El) / dE for the breakup of into a 
proton and ^Be into a s wave continuum state as a function of the cm. energy E for 
various depths Vq of the potential in the contiuum. The proton separation energy 
in MeV is given in parenthesis. 

As examples we consider the three unstable nuclei ®B [14,15,62,63,64,65,66], 
[67,68,69] and ^^N [70,71,72]. Their ground state is well described by a proton 
in a p-wave bound state. The proton wave function is calculated in a single- 
particle model. The depth of the potentials was adjusted to reproduce the 
experimental binding energies as in the case of the neutron+corc nuclei. The 
corresponding parameters of the Woods-Saxon potential are given in table 7. 
As a consequence the ANC of the bound state wave functions is uniquely 
determined in this model. 

The reduced transition probability dB{El) / dE as a function of the relative 
energy E is shown in figures 21, 22, and 23 for these nuclei for various depths 
of the continuum potential. In the calculation it was assumed that the chan- 
nel spin S of the ground state is given by 5" = jc -|- 1/2 with the core spin 
jc- The additional Coulomb interaction between nucleon and core leads to a 
substantial modification of the shape functions. The transition strength still 
peaks at low relative energies, however, the maximum is shifted to higher en- 
ergies as compared to the neutron-|-core case. For a pure Coulomb wave in 
the continuum state, i.e. Vq = 0, the maximum is reached at an energy of 
E/Sp = 4.28, 2.57 and 1.26, for ^B, ^^N, and ^C, respectively. The larger shift 
of the maximum to higher energies corresponds to an increase of the parameter 
rji, cf. Table 7. In the neutron -|-core case the maximum is expected at a value 
of only 0.18. The overall shape of the transition strength changes similarly 



48 



'N (0.601 MeV) 




Fig. 23. Same as Fig. 21 but for the breakup of into a proton and ^B. 

with the depth of the potential as for neutron+core nuclei and the occurence 
of resonances in the p-wave continuum is observed again. Also the integrated 

reduced transition probability, as depicted in figure 24, resembles in its depen- 
dence on the potential depth the results for neutron+core systems between 
the two cases of s — > p and d ^ p transitions. The B{E1) value assuming the 
same potential in the continuum as in the bound state is clearly smaller than 
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Fig. 24. Reduced transition probability B{E1) for the breakup into a proton+core 
with cm. energies E between and 10 MeV as a function of the potential depth 
Vq for the continuum states. The fiUed circle gives the result for the same potential 
depth as for the single particle bound state. The proton separation energy for each 
system is given in parenthesis. 

the value for pure Coulomb waves in the continuum. 



6 Low-energy behaviour and ANC method 



At small relative energies the reduced transition probability dB{EX)/dE oc 
Slj (A) is strongly suppressed due to penetration effects of the centrifugal bar- 
rier. It is useful to change to a different quantity to study the low-energy 
behaviour and the effects of the nucleon-core interaction in the continuum. 
Considering the approximation of the shape functions for small x in table 4 
one can use 

5;f(A)=x-^'/-^5;f(A) (76) 



with a finite limit for a; — > in the case of neutron-|-core systems. This quan- 
titiy takes the angular momentum barrier into account. Correspondingly the 
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quantities 



2lf-l 



(77) 



and 



5;/(capt, A) = x-2^'+^5;/(capt, A) = x-'^^-\l + x' f'+'Sl;{X) (78) 



are the relevant forms with a finite value for x ^ for photon absorption and 
radiative capture, respectively, cf. equations (51) and (52). (Note that for the 
capture reaction the initial and final state are interchanged.) The low-energy 
behaviour is clearly dominated by the centrifugal barrier in the continuum 
state. Prom equation (78) one finds the famous 1/v oc 1/x law of the cross 
section for the neutron capture from a s wave in the continuum Zj = 0. The 
above quantities approach a finite value in the limit E that depends on 
the function bi^.. Thus, the interaction in the continuum has an effect on the 
absolute normalization of the low-energy cross sections. 

Apart from the typical x^'^"*"^ dependence the expressions for (X) contain 
a factor (1 + x^)^^^"^ in the denominator that corresponds to a pole at x"^ = 
— 1. Considering equation (78) one sees that the shape function for capture 
reactions has a simple pole at the separation energy E = —Sn- This pole 
limits the range of convergence of an expansion of the shape functions or 
cross sections in terms of the parameter x or the energy E [46]. Since the 
separation energy is very small for halo systems a corresponding expansion 
has only limited applicability in this case. 

A similar effect is observed for proton-|-core systems [41,48]. Here, the low- 
energy behaviour is dominated by the Coulomb barrier in the continuum state. 

In nuclear astrophysics one needs, e.g., cross sections for radiative proton cap- 
ture at very low energies. Instead of extrapolating the strongly energy depen- 
dent capture cross section cr^A(p + c — > a-l-7) to low energies the astrophysical 
S factor 



is employed. It is weakly dependent on energy and approaches a finite value 
for — > 0. The exponential term depending on the Sommerfeld parameter r)i 

in the continuum state cancels in leading order the suppression of the capture 
cross section due to the Coulomb barrier. The factor E is proportional to the 
factor in the theorem of detailed balance (40) from the phase space in the 
continuum. Considering the scaling of the shape functions for proton-|-core 



S{E) = a^xip + a + ^)E exp(27r?7i) 



(79) 
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systems in section 3.3 the generalitzation of (78) is found. The quantity 

5^(capt, A) = Cl{v^)x'sl^{c^pt, A) = Cj(r;,)(l + xY'^'Sl;{X) (80) 

approaches a finite value for ,x — because Cf_{r)i) — > rjf''''^^ ex.p(—27rr]i) for 
Tji — Tjf/x ^ oo with constant rjf. It follows that 

5(E)oc5;f(capt,A)«^^ (81) 

for a; — > 0. For small separation energies of the proton in the bound state a 
strong increase of the S factor at small energies is observed due to the closeness 
of the pole aX E = —Sp. In general, the absolute value of iS'(O) depends on the 
ANC of the ground state and the scattering length, i.e. the strength of the 
continuum interaction. 

An interesting case is the direct El radiative capture from the p + ^^O con- 
tinuum to the first excited (1/2"*") state in ^^F that is relevant to nuclear 
astrophysics [44,73]. This state with I — and a small proton separation en- 
ergy of 105 keV can be considered as a typical halo state. Assuming a radius 
of R — 3.15 fm we find 7 = 0.217 and ry, = 3.787. At low energies the p ^ s 
transition gives the main contribution to the astrophysical S factor (apart 
from the small contribution of the p — > d transition to the 5/2+ ground state 
with Sp — 600 keV). Due to the closeness of the pole at £■ = —Sp a strong 
rise of the astrophysical S factor is observed for £■ — > that is practically 
independent of the potential in the continuum. 

The sensitivity of the low-energy S factor to the interaction in the continuum 
can be estimated by comparing S{0) for various p-|-core systems as a function 
of the interaction strength. In figure 25 the dependence of the S factor at zero 
energy is shown for the reactions (a) ^Be(p,7)^B, (b) ^^C(p,7)"'^^N, and (c) 
^B(p,7)^C. In all cases we have an El s ^ p transition. The S factor is cal- 
culated with varying depths of the Woods-Saxon potential in the continuum 
s waves. Spin-orbit contribution to the potential are neglected for simplicity. 
For an easier comparision all zero-energy S factors are divided by the zero- 
energy S factor for pure Coulomb waves, i.e. Vq = for the nuclear potential. 
These are given by (a) 22.07 eV b, (c) 129.4 eV b, and (b) 68.41 eV b in the 
present model. As is clearly seen in figure 25 the nuclear proton-core interac- 
tion directly has an effect on the zero-energy S factor. There is a considerable 
variation, especially if the depth Vq approaches a value where an s-wave con- 
tinuum state becomes a bound state, e.g. at Vq = 16.3 MeV and 71.6 MeV 
in case of the ^Be(p,7)^B reaction. Furthermore, we find that the sensitivity 
to the potential in the continuum state increases with the separation energy 
of the proton. Assuming the same depth for the continuum potential as for 
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Fig. 25. S factor S{E, Vq) at energy = as a function of the depth Vq of a 
Woods-Saxon potential in the continuum scaled to the S factor calculated with pure 

Coulomb waves for three radiative capture reactions. The numbers in parenthesis 
are the separation energies of the proton. The filled circle shows the result assuming 
the same potential depth in the continuum as for the bound state. 

the bound state a S factor is obtained that is somehow smaller, but close the 
value from a calculation with pure Coulomb waves. 

Our results indicate that the zero-energy S factor is not uniquely determined 
by the ground state ANC that can be extracted from experiments. In general 
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it is necessary to consider the nuclear interaction in the continuum that can be 
parametrized at small energies by the scattering lengths in the relevant partial 
waves [46,48]. One example of particular interest is the ''Be(p,7)^B reaction 
where experimental data from direct and indirect experiments for the S factor 
have to be extrapolated to zero energy assuming a certain energy dependence 
from theoretical models [74] . For the ^Be+p system s-wave scattering lengths 
as for channel spin S=2 and S=l (from coupling the spins of the proton and 
^Be) have been determined experimentally, however with large uncertainties 
[75]. More precise values are available for the mirror system '''Li+n [76]. The 
scattering lengths are easily reproduced in the single-particle model by ad- 
justing the depth of the Woods-Saxon potential in the continuum. As is seen 
clearly in figure 3 of Ref. [74] the energy dependence of the S factor shows a 
considerable variation that affects the extrapolation of the experimental data 
to zero energy. 



7 Conclusions and Outlook 

In this paper we study single-particle aspects of neutron- and proton-halo 
nuclei, i.e. the effects which appear when the separation energy of the least- 
bound nucleon tends to zero. These nuclear systems display universal features, 
e.g. for the electromagnetic transition strength, that can be described in simple 
models. 

We use different approaches in our theoretical studies: one is a square-well 
model, for which we obtain remarkably simple analytical results. A more re- 
alistic and conventional approach is the model with Woods-Saxon potentials. 
For our applications to halo nuclei we find that it is, even quantitatively, quite 
similar to the square-well model. The underlying reason is quite simple: For 
the low energies relevant for the halo systems there is shape independence: the 
properties of the potential are encoded in some low-energy parameters which 
do not depend on the details of the chosen specific potential. Halo nuclei are 
low-energy phenomena, related to large wavelengths. It is well known that a 
probe of a given wavelength is insensitive to details of structure at distances 
much smaller than this wavelength [78]. This means that we can mimic the 
real short-distance structure (e.g. determined from the nuclear many-body 
problem or a Woods-Saxon single-particle model) by a simple short distant 
structure, e.g. a square well model with its agreable analytical properties. 
The depth and radius parameters serve as adjustable parameters that, e.g., 
reproduce the binding energy and the scattering length. 

We mention also that there is a need to supplement full scale ab-initio mi- 
croscopic approaches to nuclear structure by some empirical fine tuning of 
parameters to reproduce the actual position of loosely bound single particle 
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structure, see e.g. [11]. This is indispcnsible in order to obtain a satisfactory 
description of low-lying electromagnetic strength. 

Following a study of some general single-particle properties we consider the 
radial matrix elements for electric excitations to the continuum. They essen- 
tially depend on the asymptotic wave functions of the bound and scattering 
states. We identify the important parameters and find rather simple analytical 
formulae for the most relevant cases for neutron halo nuclei. The most impor- 
tant parameter is 7 = qR depending on the separation energy S = h^q^ / (2/x) 
and the size R of the system. It is shown that 7 serves as a convenient expan- 
sion parameter. Proton halo nuclei are studied numerically. We find simple 
scaling laws for the transition strength. We also discuss typical examples for 
light halo nuclei in the Woods-Saxon model. 

We quantitatively analyze the effects of the final-state interaction in the con- 
tinuum on the transition strength. We find that it is determined mainly by the 
scattering length. In general, this nucleon-core interaction has an influence on 
the radiative capture process down to very small energies. This can serve as a 
warning for the simple application of the ANC method: the results are rather 
insensitive to the potential in the continuum only for true halo systems. 

The power of the present approach is demonstrated by recent applications to 
actual nuclei [25,26]. It will serve also as a framework for further applications, 
notably for the upcoming radioactive beam facilities RIKEN, FAIR/GSl, RIA. 
One important conclusion is that we quantitatively understand halo effects. It 
is expected from our studies that neutron-halo effects will show up irrespective 
of the mass number A, whereas proton halo effects will tend to disappear with 
increasing charge number Z. 

We hope that methods similar to ones developed in this paper will also be 
useful for the (much) more complicated problems of two-, three, and more nu- 
cleon halo nuclei [35,36]. One early discussion of that problem is [79]. A review 
on modern approaches in context of effective field theories and references are 
found in [80]. 
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Note added in proof: 



In the meantime a paper "Low-lying dipole strength for weakly bound systems: 
A simple analytic estimate" [81] appeared that contains results which are 
related to this paper. 



A Neutron in a squEire-well potential 



The problem of a neutron in a square-well potential is well studied for bound 

and scattering states in the literature, see e.g. [77], and many analytic results 
have been obtained. However, often only results for a few selected cases are pre- 
sented. Here, formulas for arbitrary values of the orbital angular momentum 
are given by generalizing earlier calculations and simplifying the expressions. 
We also derive results that, to our knowledge, were not quoted in the literature 
before. 



A.l Wave functions and probabilities 



The bound-state wave function 

Kimirl = ^YUr) (A.l) 



for a neutron with orbital angular momentum / and projection m in a square- 
well potential of Radius R and depth Vn is completely characterized by the 
separation energy S > and the principal quantum number n = 1,2,... that 
counts the number of nodes (including the node at r = 0). The radial wave 
function 



fniir) 



Anirjiiqnir) for r <R 
Bniri''hi^\iqr) for r > R 



(A.2) 



is expressed in terms of spherical Bessel and Hankel functions ji and [38] . 
The quantities q = a/2/ZS/ h and = ^2/^(1^; ~ S)/ h in the arguments are 
related by 

+ = ^ (A.3) 
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where an increasing number of radial nodes n requires larger values of qni and 
of the depth V.^^i of the potential in the ground state. The relation between q 
and for given n is determined by the continuity of the logarithmic derivative 



fni{r) 



(A.4) 



r=R 



of the radial wave function (pniir) ai r — R. This condition can be written as 
= / + 1 - = f/-) - I (A.5) 



with 



^(±) _ jl±lilnl 
— Inl 



Mini) ^'h\^\i-f) 



Z7 



(A.6) 



where 



^qniR and 7' + fni ^ ^^^2^"' ^ 



(A.7) 



The (complex) constants Ani and Bni are determined by the continuity con- 
dition for the wave function 



(A.8) 



and the normalization condition 



(A.9) 



with 



R 

\^ni? j drr"^ [ji{qnir)Y 



(A.IO) 



and 



\B.. 



nl 



dr r 



hi^\iqr) 



(A.ll) 



R 



These radial integrals are easily evaluated (see Appendix B) with the results 
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and 



(A.12) 



P> - 

^nl — 



\Bnl?' 



hf\i^)-{X,-l) 



where the continuity equation (A. 8) was used. The quantity 



X, 



r 



^(+1/2(7) 



(A.13) 



(A. 14) 



is a rational function in the variable 7. From the ratio 



P> 7^, [Xi - 1] 



(A.15) 



the probability of finding the neutron at radii r < R 

I'X, + 7^, 



Pnl — 



P< 
p< I p> 



{ill + 1') Xi 



(A.16) 



is obtained. In the limit 7 — > we have 



limX; 

7-»0 



00 for / = 
|±ifor />0 



(A.17) 



and 



lim Pn 



for i = 
l^for />0 



(A.18) 



independent of the principal quantum number n. Equivalent expressions for 
neutrons were given recently for / = 0, 1, 2 in Ref. [32]. 



In case of a scattering state the wave function has the form 



(A.19) 
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for energy E — h^k"^ / {2iJ,) with the radial wave function 



Aniknirji{knir) 



2i 



ex'p{2iSi)u\^\kr) — u\ '{kr 



(-), 



for r < R 
for r > R 



(A.20) 



where 



u\^\x) = X [-yi{x) ± iji{x)] 



(A.21) 



and ji and yi are spherical Bessel and Neumann functions, respectively [38]. 



The phase shift is denoted by 5i and kni = \j2jjL{E + Vni)/h where Vni is the 
depth of the square-well potential that gives the correct separation energy S 
of the neutron. The continuity condition for the logarithmic derivative 



^nl — 



9i{r) 



l^nlj'l{l^nl) 



= 1 + 



K 



{-)'( 



e-KY>{2i5i)ui^\K) - u\ \k) 
at r — R with 

K — kR Kfii = kjiiR 



,(-)/ 



(A.22) 



(A.23) 



determines the phase shift Si from 



exp(2i5;) 



(A.24) 



The phase shift is a sum 

5l = Ti + pi 



(A.25) 



of the hard-sphere phase shift Ti with 



exp(2iri) 



u 



(+), 



k) 



(A.26) 



and the additional phase shift pi with 



exp(2ip/) 



(A.27) 
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where 



(±)V ^ 

KUi [K) 



(A.28) 



The relation 



(A.29) 



with the depth V^i of the potential in the scattering state determines the 
quantity Rni- Generally, V^^ in (A. 7) can be different from v£. In the hmit 
k ^ we have 



tan(5„0 ^ -ai0+^ 



(A.30) 



with the scattering length 



hs (, 2/ + 1 



(A.31) 



where 



R 



(2/+l)!!(2i-l)! 



if Z = 
if Z > 



(A.32) 



is the scattering length of a hard sphere of radius R, i.e. gi{R) — correspond- 
ing to lIi^cg. 

Defining the penetrability 



Pi{x) = u^j^\x) 



X 



(A.33) 



we can write 
,(±) 



(A.34) 



with the hard-sphere phase shift ti . Then the continuity of the wave function 
at r = 
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1 r 



_i 

^Pi ^{k) exp{idi) sm{di-Ti) 



(A.35) 



fixes the constant A^i- The differential probabihty of finding the neutron inside 
the square-well potential of radius R is given by 



dPnl 
dk 



(2^ 



R 

dVtr I drr^ 



(A.36) 



TT 



7rPi(K) \ Knl Knl J 

Using the continuity relations for the logarithmic derivative one obtains 



dPni _ 7 
dx 471 



ex-p(i5i)u\^\K) — exp{—i5i)ui '(k) 



^[exp{iSi)u\fl{K) - exp{-iSi)ul_{{K)] 



(-)i 



7 / sin^ {5i -Ti) sin {Si - ti-i) sin {6i - ti+i] 



TT 



Pii"") <i ^Pi-i{k)Pi+,{k) 



(A.37) 



The probability is determined mainly by the penetrabilities P;,/±i and the 
differences of the phase shift 5i from the hard-sphere values ti ^±1 . 



A. 2 Root-mean- square radius 



The root-mean-square radius of the bound state wave function is obtained 
from 



/ 2\ -^2 



(A.38) 



with 



oo 

Rr, = j drr^+' \Mr)t 



(A.39) 
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With the help of the integrals in appendix B and the continuity equation for 
the wave function and the logarithmic derivative we find 



R 

Ro=\Anif / drr^ \ji{qnir)Y + \Bni\^ I dr h\''{iqr) 



R 



-ly 



\Bni\_ij_ 
2 



- \B 1^ 



Similarly we obtain 



2 / 1 



— + — ) Yi'-% 



(A.40) 



R 



R2^\A^ifl drr' \ji{Qnir)f + \Bj' J dr \hi'\iqr)\' (A.41) 



Qui 12 " 



R 



mm)? {Sti - '^tiY^Yt^ 



-{{2l-l)Yt'-f^,}{{2l + 3)Yr-t^}) 

-(-l)'^^[/^r^(^7)]^ (3(.7)^ - 2(^7)^>1^-VW 

- {{21 - i)i^(-) - {z^r} {{21 + 3)r/+) - {i^r}) 

i^n'i>i^^(nor|(4+^) [-2YrYr 

+ {(2Z + 3)y/+) + (2Z-i)y/-)} 

^-^){(2/-1)(2/ + 3)f/-V/+)} 



\lnl 7 

The ratio R2/R0 gives 



p2 



(-) 



(+) 



.%i T 



We note 

yj(-)yjW ^ y^(-) ^ y^(+) ^ 2Z + 1 



. (A.42) 



(A.43) 



and find 
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^ 2/ + 3 2/-1 (2/-l)(2/ + 3) 
2 H — — + - — 



1^ 



(+) 



i-2 



(+) 



ill 



(A.44) 



Explicit expressions for v/^^ are found from the recursion relations (A. 43) 
with 



(A.45) 



for Z = 0. We have 



R2 



3 + 27 ^ 3 
1 + 7 7 



/ Inl 



(A.46) 



and 



6 



5 + 57 + 272 5 5 

! ! 1 1- —— 

3 + 37 + 72 7 72;_ 



(A.47) 



for the two lowest orbital angular momenta. This leads to the divergences 
(r2)o^^ and (r^)i - ^ (A.48) 



in the limit 7^0. With the asymptotic behaviour 



Y, 



(+) 



21 + 1 for l>0 



(A.49) 



for 7 — > we obtain 



, (2/-l)(2/ + 3) , 



+ 



(2/-3)(2/ + l) 7^, 



(A.50) 



for higher orbital angular momenta. There is no divergence of the root-mean- 
square radius in the limit 7^0. It approaches a finite value. 



A. 3 Asymptotic normalization coefficient 



The modulus of the quantity B^i is obtained from the normalization condition 



1 = P<(nO + P>{nl) = d3 



\' ( 1 + ^] h\'M^l)h^'Mn) (A.51) 

Inl J 
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It is related via Bni — qCni — jCni/R to the ANC Cni- Thus we find for the 
ANC 



Cni = J2q 



■i)V i^ + i] ^l-i(n)C(^7) 

\ InlJ 



(A.52) 



In the case I — one has explicitly 

Cno = V2gexp(7) 



(A.53) 



and 



lini CnQ 

7— »o 



2q. 



(A.54) 



For Z = 1 we find 



(A.55) 



and for Z > 1 

Cni — 4/ 



2g 



(2/-3)!!(2/ + l)!! 



7' 



i-l/2 



(A.56) 



if 7 approaches zero. 



^.^ Transition integrals 



Using the notation for the radial wave functions of appendix A.l, the interior 
and exterior contribution to the radial transition integral (18) are given by 



(A.57) 



and 



l\H\>)=Bnuq 



-A-2 



exp(2i5,,)iv/;)'^(A)-A^P'^(A) 



2i 



(A.58) 
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with the dimensionless integrals 

R 



(A.59) 



and 



R 

These integrals obey the recursion relations 

+ 1 d 



<^'(A + i) = g- 

MiVi(A + l) = 
Nit'^^-\X + l)^q 



k 



dk 



If + l d 
k dk 

li + X + 2-q 



Using the continuity equations (A. 8) and (A. 35) we obtain 

il!{x,<)^A^a:i,Qnt'Mi^(^) 

^Bniiq~i 

and the ratio 



(A.60) 



Miax) 


(A.61) 


d' 

dq 


<(A) 


(A.62) 




(A.63) 


d' 
dq 


A^f^^(A). 


(A.64) 



(A.65) 



jliilnli)Knlfjlf{i^nlf 



2i 



<^(A) 



^^(A,<) 



*(A,>) 
f 7 



(A.66) 



\lnlj JlA7nk)'^nlfJl^{Knlf) [exp(2i5,^)iV;^+^'^ (A) - (A) 

With the integral formula [38] 

I dt fi{at) e gi{bt) = [afi+^{at)gi{bt) - bfi{at)gi+i{bt)] (A.67) 
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for spherical Bessel/ Neumann/ Hankel functions /; and gi we find the monopole 
functions 



M/(0) 



7 — K 



(A.68) 



and 



iVP(O) 



7^ + 



i'^'lhfU^^)uf'\K) - t'Khf\i^)u\i\{4 (A.69) 



for general orbital angular momenta I. The logarithmic derivative of the scat- 
tering wave function is given by 



(+) _ v(-) 



(A. 70) 



with 



,(-) 



jif ifinif) [exp(2i5; ^^1+^ {k) - u\j^ {k 



(A.71) 



We note 



zI-^zI+} = k' and Z^^ + Z/-^ 



21 + 1 



(A.72) 



similar as in the case for the functions i^/^^ Introducing V/^'' and Z^"^' we 
find 



K±) 



7^ — K,"^ 



(A.73) 



and 



exp(2i5,)Ar/+)'(0)-Arf^'(0) 

2 



7^ + ft;^ 



^(H-) _ ^(+) 



(A.74) 



The ratio of the monopole integrals is given by 



7^ + 



.7ni y 7ni - 7^ «nZ Vi - V j - 7^ - «;2 



(A.75) 
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where the relations (A. 7) and (A. 29) were used. If the depths of the potential 
in the bound and scattering state are identical, i.e. Vi — Vf, we have 



R'm^-l i.e. //(0) = 



(A.76) 



and the bound and scattering wave functions are orthogonal. 

Applying the recursion relations (A. 61) and (A. 63) the relevant integrals for 
El transitions / ^ / + 1 are found to be 



X 



(72 - K^f 

if - R') + zi;l) - 2r/+)z(;) + (2/ + 3)k' - (2/ + i)r 



(A.77) 



and 



exp(2^<5,+OA^^^'^'(l)-A^/'^'^'(l) 



(A.78) 



(72 + ^2)2 
X 



;(f + «^) (i^w + ztJ) + 2i^wza - {21 + iw - (21 + sy 

Then the radial integrals are 



w+i 



i?2 



;fnl{R)9nl+l{R) 



X 



(f - R') (y/+) + zl;l) - 2r/+)zfc) + {21 + 3)r' - {21 + i)f 



(A.79) 



and 



//+^(1,>) = 



X 



(72 + «;2)2 

(f + k') (y/+) + Zj-l) + 21^(+)zfcJ - (21 + 1)7^ - (21 + 3)k' 



(A.80) 



Similarly, the relevant integrals for El transitions I + 1 —>■ I 



K7"^ 



rJ/+i(7)j/(«) 



X 



(72 — fi;2)2 

{21 + 3)f - {21 + 1)R- - (f - R^) + Zt^) - 2YI^IZ^^ 



(A.81) 



and 
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exp(2i5,)A^J)'(l)-A^Vt(l) 



(A.82) 



(7' + i^^y 



(+), 



X 



{21 + 3)f + {21 + ly + 2i^y - (7^ + (y,y + 



are found with the recursion relations (A. 62) and (A. 64). Correspondingly, the 
radial integrals 



^m(l,<) 



(72 - 



;fr,i+,{R)g*jR) 



(A.83) 



X 



(2/ + 3)f - (2/ + 1)«:^ - 2F,y - (f - k^) (y^;^ + Z^^)) 



and 



/?2 



/l+i(i,>) - j:^^-^:^^fni+i{R)9*ni{R) 



X 



(-) 



(2Z + 3)7^ + {21 + + 2YiXIZl 
are obtained. 



(A.84) 



-(7^+^^) 



Assuming that the depths of the potential are the same in the bound and the 
scattering state, i.e. 



V = 'j'^ + j'^ = Vi = Vf = — 



(A.85) 



we find for the total integrals 
2R'^v 



(7' + «')' 



UR)9:i^^{R) 



(A.86) 



and 



2R^v 



(7^ + H?) 



;fni+i{R)gni{R) 



(A.87) 



consistent with the result (29) of the commutator relation. The ratios of the 
interior to the exterior dipole integral are given by 



(A.88) 
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(7^ + z^^) [yi^^^ + Z^l) + 2y/+)zti - {21 + 1)7^ - {21 + 3W 

and 



(2/ + 3)7^ - {21 + - 2y,(,',)z(+) - (7^ - R') {y^-^, + zj^^ 
{21 + 3)72 + {21 + + 2F,y - (72 + (Y^i-) + ^(+)) 

for the transitions Z — > Z + 1 and Z + 1 — > I, respectively. The scahng behaviour 
of these ratios is obtained by an expansion in terms of small 7 with k — > 0. In 
the case with a s wave bound state we have the relation 

Yi'^ = -7 = 7 cot 7 . (A.90) 



Then we find the expansion 



n n \ / n 



with Sn = (2n — l)7r/2 for 7 — > depending on the principal quantum number 
n — 1,2, Similarly, for the case of a p wave bound state the relation 

y}-^ = ^ = ^ (A.92) 

1 + 7 1 — 7 cot 7 

leads to the expansion 



7n — ^r. 



l + S-S + W-2)^ + ---) (A.93) 



with Sn — mr. Using 



— 2 

Zi~^ = and ^1-RcotR {AM) 

1 — K, cot fi; 



we obtain 



4 2 
lim = + . . . and lim = + . . . (A.95) 

n n 



for small 7. 
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B Integrals of spherical cylinder functions 

With the recursion relations [38] 



fliz) = fUz) - ^-^Mz) = -Mz) - Mz) (B.l) 
z z 

11 A- 1 

f,^,^z) + h.,{z) = ^h{z) (B.2) 

z 

for any spherical Bessel, Neumann or Hankel function fi{z) the well-known 
integral formula 

/ dz z' [fi{z)f = ^ {[fi{z)f - fUz)U,{z)) (B.3) 

is easily proven. In a similar fashion the new relation 

/ dz z' [Mz)f = ^ {3[Mz)r - 2fUz)fi+i{z) - fi-2{z)Mz)) (B.4) 

is obtained. Similarly, integrals with higher powers in z can be treated. Intro- 
ducing the logarithmic derivative 

U{z) = = / - Yl^\z) = Yt\z) -I -I (B.5) 



with 



Yl^H^) = . (B.6) 



one finds 



fi-,{z)h^,{z) = [/ + 1 + U{z)] [I - mz)] (B.7) 



and 



ii-<,{z)hmz) = {(2/ - l)[Li{z) + 1 + 1]- z'} (B.8) 

{{21 + 3)[l - Li{z)] - z'} 
With these relations we obtain 
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/ dz [Mz)f = I [Mz)f {z'-[l + l + U{z)] [I - U{z)]) (B.9) 



and 



(B.IO) 



jdzz' [Mz)f 

-{^[MzW {Sz' -2z'[l + l + L,{z)] [I - L,{z)] 

-{(21-1) [1 + 1 + Li{z)] -z''}{{2l + 3) [I - Li{z)] - z^}) 



C Radial transition integrals for neutron+core systems 

In the case when 6 is a neutron analytical expressions of the functions Ti!'/. (A) 
are found. With r] — the Coulomb wave functions 



Fiim z) zji{z) Gi{r]; z) -zyi{z) 



(C.l) 



in the integral (48) reduce to spherical Bessel and Neumann functions [38]. 
Similarly, the Whittaker function 



(2gr) ^ -qr hf\iqr) 



(C.2) 



reduces to a spherical Hankel function. This leads to 

oo 

n'/^ (A) = -i'^/.7'+2 / dt t'+' hf^ (Kt) hf^ {i^t) (C.3) 
1 

with the quantities n — kR and 7 = qR. From the recursion relation 



l__d_ 
z dz 



(C.4) 



for any spherical Bessel function gi the recursion relations 



^ti(A + l) 



If + l d 
K dn 

li + \ + 2--i 



d ' 

d'j 



(C.5) 
(C.6) 
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for the functions Ti!/. (A) are obtained that allow the calculation of all relevant 
functions for increasing values of Zj, Z/, and A from basic integrals with A = 0. 
Explicitly one finds 



dt exp[— (7 + ihi)t\ = — -j-^ exp(— 7 — m) 



(C.7) 



and 



00 / \ 

H\ (0) = / * exp[-(T- + (1 + -1 + 1 + -L, j (C.8) 



7 



K 7 + W 



exp(— 7 — in) 



with the help of the exponential integral and its recursion relations [38] . Sep- 
arating real and imaginary parts one obtains the reduced monopole integrals 



2^(0) = f [k cos{k + 5o) + 7 sin(« + <5o)] 



7^ + 

7exp(-7) 
7^ + k"^ 



/ 2 _^ ^2 \ 

-7C0s(«; -\- 5i) -\- \k-\ sin(«; + ^i) 

V 7K / 



(C.9) 
(C.IO) 



where the addition theorems of the sine and cosine functions have been used. 
The reduced integrals 2^;' (A) for fixed phase shift 5 obey the same recursion 

relations (C.5,C.6) as the functions T^/. (A). 

Considering the general form (53) for the reduced radial integrals the monopole 
functions 



(-)o. 



(C.ll) 



and 



7^^)^(o) = -K n[-^\o) = + 7) + 7I/7' 



(C.12) 



are extracted from Eqs. (C.9) and (C.IO). General expression for Tz\^^^^ (X) 
in Eq. (53) for larger values of If, k, and A can be obtained by applying the 
recursion relations 



'Rif^'-'iX + 1) = \2k'{X + 1) + (f + K'){2lf + 1)1 -Rif^iX) (C.13) 
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(A + 1) = [2f (A + 1) + (f + (i, + A + 1 - + 7) 



x7eJf'^(A)-7(f + «:^)|^7^f'^(A) 



Explicitly one finds 



(C.14) 





-«:[«:2(7-2) + 73] 


(C.15) 


ni-^\i) = 


+ «:^7(3 + 7) + 7^ 


(C.16) 


7^S+)°(l) = 


k[«:2(1 + 7)+72(3 + 7)] 


(C.17) 


M-)°(i) = 


72[k2 + 7(2 + 7)] 


(C.18) 


7^S+^'(l) = 


+ 7) + /€V(6 + 7) + 37I/7' 


(C.19) 


7^S"^'(l) = 


[/€^(3 + 37 - 7^) + ^272(6 + 7 - 72) + 37^]/72 


(C.20) 


7^w^(l) = 


-/s:[/s:2(l + 7)+72(3 + 7)] 


(C.21) 




[«;^(3 + 37 + 7^) + «;2^2^6 + 67 + 7^) + 7^(3 + 'y)]/'y^ 


(C.22) 


for A = 1 and 







7eW2(2) ^ _^[^6 ^ ^4(_g ^ 7^ ^ 27^) + k'^^^IO + 7) + 37^] 



7^S"^'(2) = fi:^(5 - 7) + A(15 + 67 - 27') + «;V(10 + 7-7') 
+37^ 

7^S+)'(2) = «;[«^(2 + 27 - 7^) + 2k'^'(5 - 7^) - 7^(2 + 7)] 



n[-^\2) = «;^(1 + 7) + «;V(7 + 27) + «:V(10 + 7 
+7^(2 + 7) 

7^(+)°(2) = n[n\3 + 37 + 7') + 2«:^=^(5 + 57 + 7') 
+7^(15 + 77 + 7^)] 

y'[K\l + 7) + 2k272(3 + 7) + 7''^(8 + 57 + 7')] 



(2) 



(C.23) 
(C.24) 

(C.25) 
(C.26) 

(C.27) 

(C.28) 



for A = 2. They cover most of the relevant El and E2 transitions between s, 
p, and d waves in the initial and final states. 
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